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FOREWORD 

This  document  is  Volume  XVII  of  the  Interim  Report  series  for  the  Pas- 
sive Nosetip  Technology  (PANT)  program.  A summary  of  the  documents  in  this 
series  prepared  to  date  is  as  follows: 

Volume  I - Program  Overview  (U) 

Volume  II  - Environment  and  Material  Response  Procedures  for  Nosetip 
Design  (U) 

Volume  III  - Surface  Roughness  Data 

Part  I - Experimental  Data 

Part  II  - Roughness  Augmented  Heating  Data  Correlation 
and  Analysis  (U) 

Part  III  - Boundary  Layer  Transition  Data  Correlation  and 
Analysis  (U) 

Volume  IV  - Heat  Transfer  and  Pressure  Distributions  on  Ablated  Shapes 
Part  I - Experimental  Data 
Part  II  - Data  Correlation 

Volume  V - Definition  of  Shape  Change  Phenomenology  from  Low  Tempera- 
ture Ablator  Experiments 

Part  I - Experimental  Data,  Series  C (Preliminary  Test 
Series) 

Part  II  - Experimental  Data,  Series  D (Final  Test  Series) 

Part  I'.I  - Shape  Change  Data  Correlation  and  Analysis 

Volume  VI  - Graphite  Ablation  Data  Correlation  and  Analysis  (U) 

Volume  VII  - Computer  tier's  Manual,  Steady-State  Analysis  of  Ablating 
Nosetips  (SAANT)  Program 

Volume  VIII  - Computer  User's  Manual,  Passive  Graphite  Ablating  Nosetip 
(PAGAN)  Program 

Volume  IX  - Unsteady  Flow  on  Ablated  Nosetip  Shapes  - PANT  Series  G 
Test  and  Analysis  Report 
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ABSTRACT 

A computer  program  is  developed  to  numerically  model  the  in-depth 
transient  response  and  shape  history  of  an  ablating  nosetip  subjected  to  a 
reentry  environment.  The  generality  of  the  input  also  allows  the  user  to 
conveniently  analyze  the  boundary  layer,  shape  change  and  in-depth  response 
of  many  materials  in  ■»  variety  of  test  facilities.  The  computer  code  is 
capable  of  handling  nosetips  of  shell  or  plug  geometries.  The  boundary  layer 
and  heat  transfer  distribution  are  modeled  for  a variety  of  environments 
including  hydrometer  erosion.  In  addition,  inviscid  flow  and  heat  transfer 
distributions  for  many  types  of  blunt  bodies  in  hypersonic  flow  can  be  readily 
calculated. 
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SECTION  1 


INTRODUCTION 

The  purpose  of  this  document  is  to  provide  a description  of  the  modeling 
techniques  and  input  requirements  of  the  EROsion  Shape  (EROS)  computer  code 
that  combines  environment  modeling  techniques  developed  by  Aerotherm  primarily 
under  the  PANT  program  (Reference  29)  with  in-depth  transient  conduction  rou- 
tines developed  at  the  Aerospace  Corporation  (Reference  28) . 

The  primary  purpose  of  this  code  is  to  numerically  model  the  in-depth 
transient  response  and  shape  history  of  an  ablating  nosetip  subject  to  a reentry 
environment.  The  code  calculates  the  inviscid  flow  and  heat  transfer  distri- 
bution for  many  types  of  blunt  bodies  in  hypersonic  flow.  In  addition,  the 
boundary  layer  and  heat  transfer  distributions  are  modeled  for  a variety  of 
environments  including  hydrometer  erosion.  The  in-depth  thermal  response  is 
capable  of  calculating  the  three-dimensional  temperature  field  and  surface 
recession  of  nosetips  at  angle  of  attack.  However,  due  to  limitations  of  the 
environment  package  to  axi-symmetric  geometries,  the  present  code  is  restricted 
to  nosetips  at  zero  angle  of  attack.  A general  thermochemistry  model,  including 
kinetic  effects,  is  used  in  the  surface  energy  balance  formalation. 

The  generality  of  the  input  allows  the  user  to  conveniently  analyze  the 
boundary  layer,  shape  change  and  in-depth  response  of  many  materials  in  a 
variety  of  test  facilities,  including  wind  tunnel,  ballistic  range,  and  arc 
heater. 

A description  of  the  numerical  modeling  and  calculation  procedure  is 
given  in  Section  2.'  Input  requirements  and  output  are  described  in  Section  3 
and  a sample  problem  is  presented  in  Section  4. 
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SECTION  2 

NUMERICAL  MODELING  AND  COMPUTATIONAL  PROCEDURES 

The  problem  modeled  by  the  computer  code  is  that  of  determining  the 
instantaneous  shape  of  an  ablating  axisymmetric  nosetip  reentering  the  atmos- 
phere at  zero  degrees  angle  of  attack,  as  well  as  the  in-depth  heat  transfer 
and  temperature  rise.  The  requirement  that  the  flow  be  parallel  to  the  body 
centerline  reduces  the  problem  to  one  of  axisymmetric  flow  and  recession.  The 
nosetip  shape  change  events  are  modeled  using  the  cyclic  calculation  procedure 
outlined  in  Figure  2-1. 

The  five  computation  elements  illustrated  in  Figure  2-1  are  described 
in  the  following  sections.  Section  2.1  covers  the  inviscid  flow  solutions; 
Section  2.2  described  the  boundary  layer  heat  and  mass  transfer  calculations; 
Section  2.3  describes  the  details  of  the  in-depth  conduction  calculation; 

Section  2.4  explains  the  surface  ablation  calculations;  and  Section  2.5  de- 
scribes the  body  movement. 

The  computational  scheme  is  stable  and  accurate  only  if  computational 
time  steps  are  kept  within  certain  limits.  These  limits  are  imposed  by  in- 
depth  conduction  and  snape  change  and  are  described  in  Section  2.3  and  2.5, 
respectively. 

2.1  INVISCID  FLOW  FIELD 

The  inviscid  flow  field  serves  as  a boundary  condition  for  the  boundary 
layer  solution.  The  actual  boundary  layer  edge  state  is  determined  from  the 
shock  shape  and  the  pressure  distribution.  The  following  three  sections  de- 
scribe the  shock  shape,  pressure  distribution,  and  boundary  layer  edge  state 
calculations,  respectively.  A complete  description  and  justification  of  the 
inviscid  flow  field  calculation  technique  is  given  in  Reference  1. 

2.1.1  Shock  Shape  Calculation 

The  bow  shock  geometry  ahead  of  an  axisymmetric  ablated  shape  is  computed 
by  forming  a piecewise  1 ' ;ar  curve  with  line  segment  slopes  and  lengths  depen- 
dent on  body  point  slopes  and  spacing.  The  technique  for  evaluating  respective 
shock  point  locations  is  given  in  Figure  2-2.  The  procedure  is  to  step  along 
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Figure  2-2.  Shock  shape  evaluation  technique. 
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the  shock  by  computing  the  next  point  based  on  the  previous  shock  point  and  the 
surface  slope  quantities.  The  following  are  needed  to  perform  the  calculation: 

• Body  geometry  quantities 

• Shock  standoff  distance  on  the  stagnation  line 

• An  expression  which  relates  local  surface  slope  to  shock  slope. 

A correlation  based  on  the  results  of  Reference  30  is  used  to  compute  the 
standoff  distance  (AQ) . The  correlation,  which  includes  dependencies  on  stream 
Mach  number,  specific  heat  ratio  and  body  bluntness  ratio  (i.e.,  t*/z*)  is 
given  by 

/A_  \ 


m 


CR  r* 
sphere 


(2-1) 


where 


ft) 


sphere 


(Y-l)  ti*  ♦ 2 + x( 
4 (Mi-1)  \ 


1 

2K  (z ) 


r*  is  the  sonic  point  ordinate  and  CR  is  a bluntness  ratio  correction  factor 
which  is  a function  of  r*/z*. 

For  subsonic  flow  in  the  shock  layer  the  relationship  between  shock  angle 
and  local  body  angle  is  obtained  from  Reference  3 and  is  given  by 


30*  + + 


(2-2) 


where  0^  and  0 ^ are  in  degrees.  The  accuracy  of  this  approach  is  discussed 
in  Reference  1. 

For  supersonic  flow  in  the  shock  layer,  the  shock  angle  is  determined  by 
the  tangent  cone  approximation.  In  this  approach  the  shock  angle  is  a function 
of  stream  Mach  number,  ratio  of  specific  heats,  and  body  angle. 

The  procedure  is  used  with  all  environment  options  except  that  for  the 
arc  heater.  Since  hosetip  shape  change  tests  in  arc  heater  environments  are 
generally  at  a relatively  low  free  stream  Mach  number,  it  is  more  accurate  to 
use  a normal  shock  assumption  for  all  inviscid  flow  calculations. 

2«1.2  Pressure  Distribution 

The  pressure  distribution  calculation  is  based  on  regional  correlations 
as  indicated  in  Figure  2-3.  Region  I is  defined  as  the  subsonic  portion  of  the 
flow  forward  of  the  limiting  sonic  flow  characteristic.  Region  II  is  the  super- 
sonic forebody.  Region  III  is  the  flow  over  the  aft  conic  surface  of  the  nose- 
tip  and  starts  at  the  point  where  the  body  slope  approaches  the  cone  half  angle. 
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A flow  chart  identifying  the  various  aspects  of  the  pressure  distribution  cal 
culation  is  given  in  Figure  2-4.  The  procedures  used  to  compute  the  pressure 
in  the  three  regions  are  described  in  the  following  subsections. 
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2. 1.2.1  Region  I,  Stagnation  Point  to  Sonic  Point 

The  correlation  of  Reference  4 as  improved  in  Reference  5 is  used  in 
this  region  to  more  realistically  represent  the  stagnation  point  velocity  gra- 
dient and  subsonic  region  pressure  distribution  on  very  blunt  bodies.  In  addi- 
tion, the  correlation  was  extended  to  include  low  free  stream  Mach  numbers. 

The  correlation  is  an  empirical  extension  and  syn thesis  of  the  modified 
Newtonian  correlation,  valid  for  spheres,  but  including  a correlation  for  flat 
faced  cylinders.  It  is  expressed  as  follows: 
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where 


P * P/P r 


Po  « stagnation  point  pressure 


*N  « stagnation  point  radius  of  curvature 
R^  = max  (1^,  R*) 

R*  = distance  from  sonic  point  to  body  axis,  measured  normal  to  the 
surface  at  the  sonic  point 

s * surface  wetted  length  from  stagnation  point 

6 » angle  local  tangent  makes  with  body  axis 

* * sonic  point 


P-  “ P«/PG 

Pmn  - P-  + (1  ~ P Jsin*9 
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(2-4) 

(2-5) 


(2-6) 


as  it  a as 


i 


ppu  - 1 - «"n(l  - P*)  - ^ (S/a*)2  - e"n  (2-7) 

(flat  faced  cylinder  pressure  distribution) 

with 

n “ 5 ^ln(sVs)  (2-8) 

The  sonic  point  location  is  an  important  parameter  in  calculating  the 
surface  pressure  distribution,  for  it  determines  the  surface  length  over  which 
the  Region  I correlation  is  used  in  the  "subsonic"  nose  region.  The  importance 
of  the  sonic  point  reflects  the  fact  that  geometry  effects  downstream  of  the 
3onic  point  have  no  influence  on  the  subsonic  region  flow  field,  and  hence, 
pressure  distribution.  In  the  code  the  sonic  point  is  found  using  correlations 
which  account  for  the  effects  of  the  following: 

• Free  stream  Mach  number  (M  ) 

00 

• Ratio  of  specific  heats  (y) 

• Nose  tip  bluntness  (r*/z*) 

• Surface  streamline  recompression  on  biconic  shapes  (B*) 

The  procedure  is  to  estimate  a sonic  point  location  assuming  modified  Newtonian 
flow  and  then  correct  the  location  for  the  effects  noted  above. 

The  modified  Newtonian  sonic  point  is  the  first  point  downstream  of  the 
stagnation  point  (6  = 0),  whith  has  an  angle  (B)  greater  than  the  following: 


where  p^  and  p*  are  defined  in  Equations  (2-4)  and  (2-6),  respectively,  and  B* 
is  defined  in  Figure  2-3. 

The  nosetip  geometry  is  then  interrogated  to  determine  the  following: 

• Bluntness  ratio  at  Newtonian  sonic  point 

• The  existence  of  a conic  surface  with  30°  <8  < 60"  and  the  conic 

c 

surface  half  angle  B 

c 

The  bluntness  ratio,  specific  heat  ratio,  and  free  stream  Mach  number  are  used 
to  obtain  a blunt  body  sonic  point  from  correlations  of  exact  numerical  pre- 
dictions . 
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M > 7 
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8*  = 51.3 
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The  < urves  are  hyperbolas.  The  expressions  are  written  in  a form  to  illustrate 
as  clearly  as  possible  their  interrelationships. 

If  a conic  surface  is  formed  in  the  vicinity  of  the  sonic  point,  the 
minimum  cone  half  angle  for  supersonic  cone  flow  is  also  computed  using  the 
following  correlations  of  exact  solutions. 
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8*  « 19.1  + 28.2*e 


..  -8.52%i 


(M  - 2.0) 


For  other  values  of  y»  linear  interpolation  is  used. 


If  the  cone  half  angle  ({$  ) is  greater  than  then  the  cone  is  super- 

C C 


sonic  and  the  sonic  point  i3  at  the  forward  end. 

The  logic  to  decide  whether  the  sonic  point  i«  controlled  by  cone  flow 
or  blunt  body  flow  is  shown  in  Figure  2-5. 


2. 1.2. 2 Kegion  II,  Sonic  Point  to  Match  Point  with  Aft  Body  Correlations 

In  t!»e  supersonic  forebody  of  the  nosetip  (Region  II)  pressure  distribu- 


tions are  computed  either  using  the  modified  Newtonian  expression  (Equation 
(2-5))  or,  for  biconic  type  configurations,  using  a conic  surface  recompression 
correlation.  The  cone  recompression  model  is  based  on  sphere/cone  and 
ellipsoid/cone  exact  solutions  performed  at  various  Mach  numbers.  The  stream- 
wise  length  required  to  obtain  the  recompression  is  given  by: 


In 


«)  - 


805(8  - 6 ) 2 - 0.22 
o 


(2-12) 


where 


6 «=  1.047  radians 

o 


sR  *=  stream  length  from  stagnation  point  to  the  end  of  cone  recompres- 
sion 


*N  = geometric  stagnation  point  radius  of  curvature 


Along  a conic  surface  starting  at  s^,  the  recorapress ion  pressures  are  given  by 
a linear  function  of  stream  distance;  i.e.. 


s **  s • 

P " Pi  + 5T^-S-  <pc  - Pi>  *>r  sL  < s < sR 


(2-13) 
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for  s > s. 
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Figure  2-5.  Flow  chart  of  logic  to  determine  nosetlp  sonic  point 
location  (subroutine  RUNLP). 
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where 

Pc  * PC/PQ  “ sharp  cone  pressure  ratio 

p^  “ P/PQ  = Pressure  ratio  at  start  of  cone,  s^ 


The  pressure  distribution  computation  in  Region  II  of  the  nosetip  must 
also  blend  together  the  results  from  the  several  correlations,  including  the 
following : 

e Region  I subsonic  flow  correlation 

e Region  II  conic  surface  reconpression,  if  any 

e Region  III  CQ  correlation,  Prandtl-Meyer  flow  or  modified  Newtonian 
(see  Section  2. 1.2.3) 

The  smoothing  is  performed  using  a weighted  average  between  an  incremen- 
tal modified  Newtonian  expression  and  a linear  aecay  expression.  For  smoothing 
in  the  region  9initial  > 0 > 0£inal 
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incremental 

modified 

Newtonian 


linear  interpolation 
between  end  of  merging 
region  and  sonic  point 


(2-14) 


where  a is  the  weighting  function.  Taking  a linear  weighting. 


gives 
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P + gj - (1  - p„)  (sin*0  - sinJ0i) 


(Pf  - Pi> 


(2-15) 


(2-16) 
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In  a typical  case,  the  smoothing  expression  might  be  used  between  the 
sonic  point  and  the  start  of  a forecone  surface  and  between  the  end  of  the 
forecone  and  the  match  point  with  the  Region  III  correlations.  In  the  case 
where  concave  shapes  develop  as  in  the  sketch  below,  smoothing  is  performed 
between  the  sonic  point  and  the  inflection  point  in  the  shape.  Downstream  of 
the  inflection  point,  the  modified  Newtonian  relation  (Equation  (2-5))  is  used. 

2. 1.2. 3 Region  III  - Aft  Body 

The  correlation  for  aft  cone  pressures  is  one  developed  at  Aerospace 
Corporation  (Reference  6) . It  has  the  form 

s f fls  ts. L_ 

V n\D  & 

where 

0c  = cone  half  angle 

z axial  distance  from  the  start  of 
c 

D = diameter  at  start  of  aft  cone 
CD  = drag  coefficient  of  the  forebody 

cp  " <P  ' P.)/(1/2)PX 

The  function  fR  is  determined  by  a series  of 
numerical  solution  for  cones  of  varying  bluntness,  with  cone  half  angle  as  a 
parameter.  The  curves  asymptotically  approach  the  sharp  cone  pressure  level. 

Thi'  transition  between  Regions  II  and  III  is  effected  at  the  point  where 
the  pressure  distribution  curves  for  the  two  regions  intersect.  That  point  is 
determined  'iteratively  since  CQ  is  a function  of  its  location. 

The  calculation  procedures  used  in  Region  III  (aft  of  shoulder)  are 
based  on  hypersonic  considerations.  They  are  used  for  Mw  5.  To  better  model 
the  flow  for  < 4,  the  modified  Newtonian  calculation  procedure  for  Region  II 
is  extended  to  Region  III.  For  4 < < 5 , linear  interpolation  is  used  between 

= 4 and  **  5 predictions. 

Alternate  procedures  are  used  for  cylindrical  afterbodies. 

2.2.3  Boundary  Layer  Edge  State 

The  actual  boundary  layer  edge  thermodynamic  state  is  determined  by  a 
look-up  on  pressure  and  entropy  in  a real  gas  Nollier  air  table.  Pressure  is 
known  from  the  inviscid  flow  solution,  and  entropy  is  calculated  from  considerations 
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the  aft  cone 


polynomial  curve  fits  of  exact 
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of  bow  shock  shape  and  boundary  layer  mass  flux.  The  sketch  shown  on  the 
following  page  illustrates  the  path  of  a streamline  passing  through  the  shock 
layer.  At  the  point  where  the  streamline,  originating  at  y,  enters  the  bound- 
ary layer,  the  mass  flux  can  be  expressed  as  follows: 


Boundary  layer  edge 


L 


For  the  laminar  boundary  layer 


P ■ 4*52  rPeRe0, 


(2-18) 


For  the  composite  model  of  the  turbulent  boundary  layer,  which  is  described  in 
Section  2.2.2,  the  free  stream  location  of  the  boundary  layer  edge  streamline 
is  computed  from 


p.u«y£ 


(100  + 2Re0  V 
100  + ^e) 


4.52nieRe 


(2-19) 


The  turbulent  Reynolds  number,  Reg^,  is  computed  using  a roughness  augmented 
momentum  thickness.  This  expression  for  y?  passes  smoothly  from  the  laminar 
value  at  Re0  = 0 to  the  previously  used  expression  for  turbulent  flow  at  large 

The  entropy  used  to  compute  the  edge  conditions  when  the  streamline 
enters  the  boundary  layer  is  the  entropy  existing  at  the  radial  coordinate  yL 
or  yT  just  behind  the  shock.  This  entropy  is  evaluated  from  the  free  stream 
conditions,  using  oblique  shock  relations  (see  Reference  6). 

The  boundary  layer  edge  velocity  over  most  of  the  body  is  computed  from 
energy  conservation  along  an  effective  inviscid  flow  stream  tube  as  follows. 

u - J2(H  - h ) (2-20) 
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Since  the  boundary  layer  edge  state  is  determined  accounting  for  entropy 
swallowing,  the  edge  velocity  also  is  affected.  This  is  the  only  mechanism 
by  which  swallowing  influences  the  boundary  layer  solution. 

In  the  vicinity  of  the  stagnation  point  the  velocity  is  assumed  to  be 
linear  for  Pe/Pt2  = 1,0  to  0*999  • 1116  velocity  at  the  0.999  point  is  calcu- 

lated accounting  for  entropy  swallowing  using  a first  guess  obtained  by  assum- 
ing normal  3hock  entropy.  Therefore,  the  velocity  gradient,  due/ds|  , is  eval- 
uated directly  from  the  pressure  distribution  including  entropy  layer  effects. 

The  edge  viscosity  is  determined  from  the  following  correlation  taken 
from  Reference  6. 
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T < 2000°R 


T > 2000°R 
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T in  °R,  y in  lb/ft-sec  units 
The  Prandtl  number  is  assumed  constant  at  0.7. 

2.2  BOUNDARY  LAYER  HEAT  AND  MASS  TRANSPORT 

The  boundary  layer  heat  and  mass  transport  events  are  modeled  using  a 
film  coefficient  approach.  The  momentum  integral  equation  is  solved  assuming 
that  zero  pressure  gradient  relations  between  skin  friction  and  momentum  thick- 
ness apply  in  the  presence  of  pressure  gradients.  Reynolds  analogy  and  compres- 
sibility corrections  are  applied  to  obtain  the  nonblown  heat  transfer  coefficient 
distribution.  Effects  of  blowing  are  accounted  for  as  a function  of  local  abla- 
tion rate,  and  the  mass  transfer  coefficient  is  taken  as  a constant  ratio  of 
heat  coefficients. 

Details  of  the  solution  procedure  for  laminar  flow  are  given  in  Section 
2.2.1.  The  turbulent  solution  procedure  is  discussed  in  Section  2.2.2.  Tran- 
sition criteria  options  in  the  code  are  given  in  Section  2.2.3.  Techniques 
used  to  compute  the  roughness  effects  on  laminar  and  turbulent  heating  are 
reviewed  in  Section  2.2.4,  and  relations  used  to  compute  heat  transfer  in  re- 
gions of  transitional  boundary  layer  flow  are  discussed  in  Section  2.2.5.  The 
effects  of  hydrometer  boundary  layer  stirring  are  covered  in  Section  2.2.6. 
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2.2.1  Laminar  Heat  Transfer,  Smooth  Wall 

The  stagnation  point  heat  transfer  coefficient  calculation  is  discussed 
in  Section  2.2.1. 1,  and  the  laminar  distribution  evaluation  technique  is  de- 
scribed in  Section  2. 2.1.2. 

2. 2.1.1  Stagnation  Point  Heat  Transfer  Coefficient 

At  high  altitude  or  low  Reynolds  number  conditions,  the  energy  flux  to 
the  surface  is  limited  by  the  total  energy  content  of  the  free  stream.  The 
corresponding  heat  transfer  coefficient  is,  therefore. 


qlimit  " P-U-Hc 


(2-22) 


p u C„  - = p u 

6 e H'°  limit  " " 


For  other  conditions,  the  stagnation  point  heat  transfer  coefficient  is 
computed  using  the  relation  of  Fay  and  Riddell  (Reference  7)  with  Pr  **  0.7. 


p u C„  _ * 0 
e H,0 


where,  as  suggested  in  Reference  8, 


i-»  1-01  ^ 12-“! 


!0  , T©  < 5000“R  I 

1 - 0.308(Tq/Ho) , TQ  > 5000°R  ) 


(2-24) 


The  Lewis  number  used  is  the  average  between  the  Lewis  numbers  evaluated  at 
the  wall  and  edge  temperatures.  These  are  computed  from  the  approximation 


(1.2  ■ , 

(1.2  - 5.5  x 10_*(T  - 5400°R)  , 


T < 5400*R 
T > 5400*R 


(2-25) 


2. 2. 1.2  Laminar  Heating  Distribution 

The  method  of  Reference  9 as  simplified  in  Reference  10  is  used  to  obtain 
the  laminar  heating  distribution.  The  correlation  is  expressed  as  the  local 
heat  transfer  coefficient  divided  by  the  stagnation  coefficient,  i.e.. 
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The  corresponding  laminar  momentum  thickness  Reynolds  number  is  obtained  by 
applying  Lees  transformation  to  the  Blasius  incompressible  flat-plate  skin  fric- 
tion relation.  The  resulting  equation  is 


_ 0.664 
R<\  ' 


Si  r^SL  hbi)  ur^dsT  1 
r Lpe/po l 'p°l e 


(2-28) 


This  Reynolds  number  and  the  associated  momentum  thickness  are  used  to  obtain 
boundary  layer  thickness  parameters  for  use  with  transition  criteria,  transi- 
tional heating  correlations,  and  turbulent  boundary  layer  starting  conditions. 

2.2.2  Turbulent  Heat  Transfer,  Smooth  Wall 

The  compressible  boundary  layer  momentum  integral  equation  is  solved  to 
evaluate  the  fully  turbulent  heat  transfer  coefficient  distribution.  The  im- 
portant assumptions  included  in  the  solution  are  as  follows: 

• Blowing  effects  may  be  decoupled  from  the  boundary  layer  solution 
(computes  nonblown  transfer  coefficient). 

• Boundary  layer  shape  factor  H * 6*/0  is  taken  as  -1. 
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• A modified  Reynolds  analogy  (explained  more  completely  below)  is 
used  to  relate  heat  transfer  to  skin  friction. 

• The  Crowell  incompressible  composite  skin  friction  expression  (Ref- 
erence 11)  modified  for  compressibility,  is  used. 

The  integral  momentum  equation  may  be  written  as 


S-  u,  U 6)  = i - p U e » * is  ♦ I “E 

M we  e u Me  e u ds  r ds 

e e 


(2-29) 


Using  properties  (p,  y)  evaluated  at  the  Eckert  reference  enthalpy  (Reference  12) 


H * 0.5  h„  + 0.3  h + 0.2  H 
w e o 


(2-30) 


The  Crowell  composite  skin  friction  expression  modified  for  compressibility  is 


where 
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100  + Re 
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This  expression  substituted  into  the  integral  momentum  equation  (which  is  then 
integrated)  determines  9{s).  Trial  calculations  for  sphere  cone  geometries  were 
carried  out  assuming  -1.0  < H < 0.5.  Although  H « 0.5  is  probably  most  realistic 
for  conditions  of  interest,  the  skin  friction  was  found  to  be  relatively  insen- 
sitive to  variations  in  K.  The  closest  approximation  (wxthin  10  percent)  be- 
tween 0(0)  from  the  composite  model  and  the  laminar  calculations  is  obtained  for 
H **  -1.  This  value.,  often  assumed  because  it  simplifies  the  momentum  equation, 
was  adopted  here. 

In  using  Reynolds  analogy  to  determine  the  Stanton  number  from  the  skin 
friction,  separate  factors  are  used  to  multiply  the  laminar  and  turbulent  con- 
tributions to  the  skin  friction.  The  laminar  Reynolds  analogy  factor  (taken  to 
be  independent  of  s)  is  determined  from  the  requirement  that  the  composite  model 
yields  the  correct  heat  transfer  at  the  stagnation  point,  that  is. 


n,J.*«S!jawvWW  .)>..■  itwa  a»rM.  v.sg>y..?,' j.J  rrr 


\ 


te) 


lam^  8*0 


’•J78  (t) 


(2-32) 


The  turbulent  Reynolds  analogy  factor,  F^,,  is  currently  taken  to  be  unity  over 
the  entire  body.  There  is,  however,  some  evidence  from  turbulent  BLIMP  solutions 
that  RT  is  a function  of  pressure  gradient,  being  about  0.95  on  the  nose  and 
about  1.15  on  the  cone,  therefore,  the  factor,  R^ , was  chosen  to  be 

RT  = 1 (2-33) 


This  is  close  to  the  suggestion  of  Kays'  (Reference  14)  in  which  the  exponent  on 
the  Prandtl  number  is  -0.4.  The  con^osite  expression  for  the  turbulent  heat 
transfer  coefficient  becomes 
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(2-34) 


The  turbulent  heating  factor  is  defined  as 


p « 

T peUeCH,0 


(2-35) 


2.2.3  Transition  Criteria 

Built  into  the  code  are  several  optional  techniques  for  determining  the 
conditions  for  boundary  layer  transition.  In  summary  these  are: 

• Momentum  thickness  Reynolds  number  versus  boundary  layer  edge  Mach 
number. 

• Run  length  Reynolds  number  versus  boundary  layer  edge  Mach  number. 

• Axial  distance  from  stagnation  po'nt  versus  vehicle  altitude. 

• Roughwall  transition  criterion  based  on  momentum  thickness  and  wall 
cooling  ratio. 

• Roughwall  transition  criterion  based  on  displacement  thickness  and 
run  length. 

• Fully  turbulent  flow  from  the  stagnation  point  (composite  model) . 

The  appropriate  boundary  layer  quantities  are  computed  and  compared  to  critical 
values  input  in  tabular  form  by  the  user.  For  the  two  roughwall  criteria. 
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critical  values  of  the  appropriate  parameters  from  Reference  15  jure  built  into 
the  code.  For  the  transition  correlating  parameter  involving  the  momentum 
thicknesc,  critical  values  are 
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(2-36) 


For  the  transition  correlating  parameter  involving  the  displacement  thickness, 
critical  values  are 
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where  is  the  intrinsic  roughness  of  the  surface  appropriate  for  laminar 
flow  conditions  as  specified  by  the  user  and  Re^  = Peue^i/we* 

The  onset  conditions  are  determined  first  and,  if  satisfied,  the  point 
of  transition  is  found  from  the  location  condition.  Equation  (2*3?)  indirectly 
accounts  for  the  effects  of  surface  temperature  on  transition  through  the  dis- 
placement thickness,  6*.  The  displacement  thickness  is  computed  from  the  momen- 
tum thickness,  0,  and  wall  to  edge  temperature  ratio,  (T^/T^  as  follows: 


0HOT  WALL  “ 0COLD  WALL 


[l. 104  - 0.348  (T^/T^j 


(0COLD  WALL  fr°m  Eq*  {2"28)) 


(2-38) 
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2.2.4  Surface  Roughness  Effects  on  Heat  Transfer 


Correlations  from  PANT  wind  tunnel  data  (Reference  16)  are  included  to 
account  for  the  effects  of  roughness  on  laminar  and  turbulent  heat  transfer; 
in  addition,  surface  roughness  modeling  accounting  for  intrinsic  roughness, 
scallop  roughness,  and  crater  roughness  are  included.  Roughness  effects  on 
laminar  and  turbulent  heating  are  discussed  in  Sections  2. 2.4.1  and  2.2.4. 2 
and  surface  roughness  modeling  is  covered  in  Section  2.2.4. 3. 


2. 2. 4.1  Roughness  Effects  on  Laminar  Heating 

She  effects  of  roughness  on  laminar  heating  are  correlated  with  the 
parameter 
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The  effect  is  accounted  for  with  the  multiplicative  factor  or.  the  smooth 
wall  laminar  Stanton  number,  Cu  T , where 
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The  correlation  is  applied  to  all  laminar  flow  locations. 

2. 2. 4. 2 Roughness  Effects  on  Turbulent  Heating 

The  effects  of  roughness  on  turbulent  heating  are  corr_lated  in  Refer- 
ence 16  using  the  following  parameter: 
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As  with  laminar  flow,  the  effect  of  roughness  on  turbulent  heating  is  accounted 
for  with  a multiplicative  factor  on  the  smooth  wall  Stanton  number,  CH  T« 

The  correlation  equation  is  as  follows: 
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In  the  expression  for  the  composite  turbulcrt  heat  transfer  cofficient,  the 
laminar  and  turbulent  contributions  are  augmented  individually,  so  that  cn  a 
rough  wall. 
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For  the  purpose  of  output  a composite  turbulent  augmentation  factor  is  defined 

as 
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2.2.4. 3 Surface  Roughness  Modeling 

Three  types  of  surface  roughness  are  modeled 

• Intrinsic  roughness,  k^ 

• Turbulent  or  scallop  roughness,  kfc 

• Crater  roughness,  kc 

Intrinsic  roughness  (k^  is  that  associated  with  the  basic  material  granularity 
and  is  input  as  a constant  for  each  material.  The  intrinsic  roughness  is  used 
in  the  transition  and  laminar  heating  correlations,  unless  there  are  particle 
impacts  as  described  below. 

The  turbulent  roughness  (kt)  is  the  effective  sand  grain  roughness  that 
results  from  turbulent  ablation  and  is  used  in  the  turbulent  heating  correla- 
tions. The  roughness  height,  k^,  is  specified  in  one  of  two  ways.  A uniform 
value  of  kfc  for  all  turbulent  regions  may  be  input  by  the  user;  or  a value  may 
be  obtained  by  using  the  scallop  dimension  correlation  from  Reference  17.  From 
the  correlation,  the  effective  turbiilent  region  scallop  depth  is  computed  as 
fol lows ; 
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where 

kfc  « the  effective  sand  grain  roughness  height  suitable  for  use  in 
Equation  (2-41) 

Kj  * a material  dependent  property  determined  from  experimental  data 
and  input  by  the  user  (a  nominal  value  for  graphite  is 

^ a * « 0.77. 

* 0.93  m-psia  ). 
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p * instantaneous  local  edge  pressure 


Crater  roughness  (kc)  results  from  the  inpact  of  hydrometer  particles. 


Presently,  the  assumption  of  hemispherical  craters  is  used  in  conjunction  with 
the  mass  loss  parameter,  G,  described  in  Section  2.4.2.  Therefore,  the  crater 
depth  is  derived  from 


IS!" 


(2-46) 


where 


Pp  * hydrometer  particle  density 


p * surface  material  density 
zn 


dp  * hydrometer  particle  diameter 


G * ine/m^n  * mass  loss  parameter 


k * crater  depth 
c 


Since  crater  roughness  (kc)  occurs  over  the  entire  body,  the  local  roughness 


(either  intrinsic  (k^  or  turbulent  (k^) ) is  compared  to  the  crater  roughen3s 


(k  ) and  the  larger  is  used, 
c 


2.2.5  Transitional  Boundary  Layer  Heat  Transfer 


Transitional  heating  is  computed  using  a modified  version  of  the  corre- 
lation of  Reference  18.  The  correlation  used  is  expressed  as  follows: 


CH, TRANS, R = CH,T,R  “ ^TR^e 


(2-47) 


The  values  of  ATR  and  n are  computed  differently  depending  on  the  approach  to 
fully  turbulent  heating;  that  is  inr 


1 > CH,T,R  ' CH, TRANS, R 
“ (CH,T,R  “ CH,L,r)tr 


n ” 0.85 


*TR  = [Ree  (CH,T,R  ' CH,L,r)], 


where  TR  denotes  the  values  at  the  point  of  transition.  For 


0 4 > ^/TtR  ~ CH , TRAMS , R 
“ (Sa,T,R  “ CH,L,r)__ 


n “ 2.0 


*TR  * [Ree  (CH,T,R  " CH, TRANS, r)]  q 4 


(2-49) 


where  the  subscript,  0.4,  denotes  the  point  where  A-^  is  reevaluated. 

Since  the  boundary  layer  is  transitional,  the  momentum  thickness  Reynolds 
number  (Re„)  in  Equations  (2-47)  and  (2-49)  is  computed  by  integrating  the  re- 

D 

duced  form  (i.e.,  H = -1)  of  the  momentum  integral  equation.  Equation  (2-29), 
assuming  unity  Reynolds  analogy  factor  and  using  the  rough  wall  Stanton  number 
(C„  _)  from  the  previous  boundary  layer  integration  station,  i.e., 
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It  should  be  noted  that  use  of  rough  wall  Stanton  numbers,  Cu  T „ and 

H f AJ  f K 

C„  _ 0,  in  the  transitional  heating  correlation  provides  for  a reasonable 
n#T#R 

transformation  from  the  laminar  to  the  turbulent  roughness  effects  models 
described  in  Section  2.2.4. 

2.2.6  Hydrometer  Boundary  Layer  Stirring  Effects 

Experiments  indicate  that  in  regions  of  laminar  flow  hydrometer  parti- 
cle impaction  and  subsequent  erosion  can  cause  significant  augmentation  to 
the  undisturbed  laminar  heat  transfer  rate.  An  option  is  provided  in  the 
code  to  model  this  laminar  stirring  augmentation.  The  correlation  is  in  terms 
of  the  ratio  of  the  disturbed  (stirred)  heat  transfer  coefficient  to  the  undis- 
turbed coefficient.  The  correlation  is  of  the  form 


STIRRED 
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(2-51) 
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where 

Pp  = particle  density 

= freestream  air  density 

G = erosion  mass  loss  parameter  (described  in  Section  2.4.2) 

0 «=  local  body  angle  (0  = 90°  at  the  stagnation  point) 

The  constants  (C  and  c)  are  presumed  to  be  a function  of  the  surface  material. 
Reference  19  indicates  that  graphite  data  are  best  correlated  by 

C = 0.098 

c = 0.317 

The  implementation  of  the  stirring  augmentation  logic  is  flagged  by  the 
J ROUGH  flag  described  in  Section  3.1.8.  When  the  stirring  augmentation  corre- 
lation is  employed  the  augmentation  factor  calculated  from  Equation  (2-51)  is 
compared  with  the  factor  calculated  from  Equation  (2-40)  and  the  larger  is  used. 

2.3  IN-DEFTH  CONDUCTION  CALCULATIONS 

This  section  briefly  describes  the  numerical  technique  used  to  solve  the 
heat  conduction  equation  inside  the  nosetip  and  the  coupling  between  the  sur- 
face energy  balance  relations  (Section  2.4)  and  the  in-depth  conduction  solu- 
tion. The  details  of  the  conduction  package  are  described  fully  by  Crowell 
(Reference  28).  In  this  section  only  a brief  review  of  the  procedure  is 
presented. 

Section  2.3.1  describes  the  coordinate  systems  and  governing  equations. 

The  finite-difference  formulations  of  the  differential  equations  and  their  solu- 
tions are  explained  in  Section  2.3.2.  The  conduction  time  step  control  is  dis- 
cussed in  Section  2.3.3. 

2.3.1  Coordinate  System,  Governing  Equations 

The  ir -depth  coordinate  systems  for  shell  and  plug  geometries  are 
illustrated  in  Figure  2-6.  For  the  shell  geometry  a body  oriented  coordinate 
system  which  is  located  on  the  internal  contour  is  used  (x,y,^) . For  the  plug, 
the  geometry  is  split  into  two  sections,  separated  by  a fictitious  boundary 
(shown  as  a dashed  line  in  the  figure) . The  location  of  the  interface  between 
the  two  sections  is  chosen  such  that  the  geometry  of  section  I is  exactly  that 
of  the  T'hell  configuration.  Thus,  the  coordinate  system  for  region  I is  also 
body  oriented  and  located  on  the  fictitious  boundary.  The  shape  of  the  interface 
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Figure  2-6.  Nosetlp  geometry. 
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is  taken  to  be  spherical  for  convenience.  Cylinderical  coordinates  are  used 
in  region  II  (the  shank  portion  of  the  plug) . 

In  the  body  oriented  coordinate  system  (region  I)  the  heat  conduction 
equation  for  temperature  dependent  properties  and  isotropic  conductivities 
may  be  written  as 
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where  A is  the  curvature  of  the  internal  contour  or  the  fictitious  interface. 

For  the  cylinderical  coordinates  (region  II)  the  conduction  equation 
takes  the  following  form 
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It  should  be  noted  that  due  to  the  axisymnetric  assumption,  there  is  •'o  tem- 
perature variation  in  the  ^-direction,  although  the  conduction  package  is 
capable  of  handling  full  three-dimensional  problems. 

The  boundary  conditions  on  all  surfaces  consist  of  specified  heat  flux. 
For  most  nosetip  problems  all  the  boundaries  except  the  receding  surface  are 
insulated  and  these  fluxes  are  zero.  In  case  of  the  plug  geometry  the  ficti- 
tious interface  between  regions  I and  II  is  not  a boundary  of  specified  flux 
or  temperature.  The  temperature  distribution  along  this  boundary  is  computed 
by  requiring  that  the  temperature  and  heat  flux  in  the  regions  I and  II  be 
identical  at  the  interface.  At  the  receding  surface  the  boundary  condition  is 
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where  n denotes  the  direction  normal  to  the  nosetip  surface  and  the  functional 
form  of  is  determined  from  the  surface  energy  balance  formulation 

(Section  2.4.1). 


In  order  to  solve  the  moving  boundary  conduction  problem  over  a fixed 
domain,  the  surface  movement  is  incorporated  into  the  heat  conduction  equation 
through  the  use  of  the  following  transformations 
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Region  I: 
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Region  II: 
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The  differential  equations  (2-525  and  (2-53)  and  the  boundary  conditions 
(2-55)  are  transformed  into  the  new  coordinate  system  and  then  solved  by  a 
finite-difference  technique.  The  description  of  this  finite-difference  pro- 
cedure is  given  in  the  following  section. 


2.3.2  Finite-Difference  Formulations 

The  finite-difference  scheme  which  has  been  adopted  is  the  Dufort- 

Frankel  method  that  is  an  unconditionally  stable  explicit  technique.  This 

method  uses  a central  time  difference  and  therefore,  requires  storage  at  two 

time  levels.  The  Dufort-Frankel  method  does  not  require  the  restrictive  lime 

Ax2 

step  limitation  of  standard  explicit  technique  (At  < -j^-) , but  for  consistency 
purposes  it  requires  that  At  goes  to  zero  faster  than  Ax. 


Variable  mesh  spacing  is  used  throughbut.  First  order  derivatives  are 
written  in  second  order  central  or  one  sided  difference  forms  and  for  the 
second  order  terras  the  Dufort-Frankel  form  is  used.  For  the  details  of  the 
finite-difference  formulation  the  reader  is  referred  to  Reference  28. 


When  the  equations  are  differenced,  the  left  hand  side  will  contain  the 

temperature  of  a node  (jkl)  at  the  n+1  time  level  and  the  right  hand  side  will 

contain  the  temperatures  of  the  neighboring  nodes  at  n-1  and  n time  levels  and 

„n+l 


The  difference  equation  is  then  solved  for  T^^* 


known  geometric  parameters. 

In  order  to  start  the  calculntions  both  the  n-1  and  n time  levels  are  set  equal 
to  the  initial  temperature  distribution. 


In  region  I,  the  domain  over  which  the  temperatures  are  obtained  from 
the  differential  equations  runs  from  J = 2,  JMAX-1  in  the  X-direction,  K = 2 to 
KMAX-1  in  the  Y-direction  and  L = 1 to  LMAX  in  the  4-direction.  In  region  II, 
the  temperatures  are  calculated  from  the  differential  equations  for  JP  «=  2 to 
JPMAX-1,  KP  = 2 to  KPMAX-1  and  L = 1 to  LMAX  in  X,  Z and  4~directions  respec- 
tively. In  the  present  axisymmetric  code,  the  computations  are  only  performed 
in  the  L = 1 plane.  The  boundary  temperatures  are  calculated  from  the  finite- 
difference  forms  of  the  boundary  conditions. 
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Following  the  calculations  of  the  interior  point  temperatures  from  the 
difference  equations,  the  centerline  values  (X«0,  J*1  and  JP»1)  are  obtained 
by  back  extrapolation  from  the  known  values  of  J=2  and  J*3  nodes.  For  axi- 
symmetric  nosetips  the  following  condition  must  be  satisfied  along  the 
centerline: 

H - 0 (2-55) 

These  derivatives  are  written  in  one-sided  forward  difference  forms  and  solved 
for  centerlir.e  temperatures  including  the  stagnation  point  temperature. 

The  surface  temperatures  and  recession  rates  are  determined  from  simul- 
taneous solution  of  tht  difference  form  of  Equation  (2-55)  with  the  surface 
energy  balance  relations. 

2.3.3  Time  Step  Control 

The  computational  time  steps  are  controlled  by  a comprehensive  technique 
to  achieve  numerical  stability,  economy  and  output  versatility.  The  code  has, 
basically,  two  kinds  of  time  steps:  a conduction  time  step  and  an  environment 

time  step.  The  print  time  step  is  currently  set  equal  to  the  environment  time 
step. 


2. 3. 3.1  Conduction  Time  Step 

The  time  step  of  conduction  calculations  is  t>-  minimum  of  the  following 

values : 

• Explicit  finite-difference  stability  limit:  d2/4a  where  d is  the 

minimum  mesh  size  and  a is  the  thermal  diffusivity  of  the  nosetip 
material.  This  is  not  currently  in  use  because  the  Dufort-Frankel 
scheme  is  unconditionally  stable. 

• Surface  temperature  rise  control:  At^^fSTRD/STRM)  where  STRD  is 

the  input  desired  surface  temperature  rise  in  one  time  step,  and 
STRM  is  the  maximum  surface  temperature  rise  achieved  during  the 
previous  time  step. 

• Surface  heat  flux  rise  control:  AtQ^d (qold/qnew)CTF,  where  *3 

the  maximum  surface  heat  flux  and  CTF  is  the  desired  growth  factor. 

A recommended  expression  for  CTF  in  terms  of  the  desired  maximum 
surface  temperature  rise  is  the  following: 

CTF  » 1.5  + (STRD  - 140) /300  (2-56) 
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• Surface  recession  control:  5/S . where  S „ is  the  maximum  value 

niax  max 

of  surface  recession  rate  and  5 is  the  smallest  distance  between  the 
first  and  second  nodes  in  the  Y-direction. 

At  the  first  conduction  step  when  a majority  of  the  above  quantities  cannot  be 
calculated,  the  following  time  step  is  also  used. 

2. 3. 4.2  Environment  Time  Step 

The  environment  time  step  determines  the  frequency  with  which  the  invis- 
cid  and  viscid  solutions  are  to  be  updated  and  is  equal  to  the  user  specified 
value  in  the  absence  of  time  step  stability  criteria.  Xn  the  presence  of  sta- 
bility criteria,  the  environment  is  redefined  whenever  either  one  or  both  of 
the  following  conditions  are  satisfied. 

• If  the  local  surface  temperature  changes  by  a factor  greater  than 
1.4. 

• If  the  tangent  of  the  local  body  angle  changes  by  a factor  of  two 
or  more. 

The  reference  condition  for  the  above  two  tests  is  the  last  environment  defini- 
tion. 

The  computation  time  step,  DTH,  is  the  minimum  of  the  conduction  and  the 
environment  time  steps.  Furthermore,  the  computation  is  automatically  termin- 
ated if  the  computed  time  step  is  less  than  the  user  specified  minimum,  DLTMIN. 

2.4  SURFACE  ABLATION  RESPONSE  CALCULATIONS 

The  formulation  of  the  surface  energy  balance  technique  used  to  compute 
the  surface  ablation  response  is  discussed  in  Section  2.4.1;  modeling  of  erosion 
due  to  hydrometer  impacts  is  described  in  Section  2.4.2;  computer  codes  to  pro- 
vide the  necessary  input  data  are  described  in  Section  2.4.3;  and  simplified 
means  of  the  surface  energy  balance  equation  are  presented  in  Section  2.4.4. 

2.4.1  Surface  Energy  Balance  Formulation 

The  ablation  rate  and  surface  temperature  at  points  on  the  nosetip  are 
determined  by  accounting  for  energy,  mass,  and  species  conservation  at  the 
ablating  surface.  The  sketch  below  illustrates  the  ablating  surface  control 
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Sketch  of  Surface  Energy  Balance  Control  Volume 
and  Energy  Flux  Terms 


volume  and  the  energy  fluxes  of  interest.  The  surface  energy  balance  equation 
employed  is  of  the  convective  tranafer  coefficient  type.  In  the  program,  this 
energy  balance  equation  takes  the  following  form: 
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(2-58) 


Before  commencing  a term  by  term  discussion  of  Equation  (2-58) , however, 
it  will  be  useful  to  describe  the  general  nature  of  this  transfer  coefficient 
expression.  Like  all  such  expression.  Equation  (2-58)  is  an  approximation,  the 
usefulness  of  which  depends  mainly  on  the  validity  of  the  transfer  coefficient 
approach.  A discussion  of  this  subject  is  far  beyond  the  scope  of  the  present 
document.  It  may  be  observed  here  that  transfer  coefficients  have  successfully 
correlated  both  data  and  "exact"  solutions  in  simple  heat  or  mass  transfer  prob- 
lems, and  in  combined  heat  and  mass  transfer  problems  for  unity  (or  near  unity) 
Lewis  number.  Equation  (2-58)  attempts  to  extend  the  transfer  coefficient 
approach  to  both  nonunity  Lewis  number  and  unequal  mass  diffusion  coefficient 
problems,  still  allowing  for  chemical  reactions  and  net  mass  transfer  effects. 
This  approach  was  suggested  in  Reference  20.  Its  validity  is  discussed  in 
References  21  and  22. 

In  Equation  (2-58) , the  term  q represents  the  "sensible  convective 

SCIl 

heat  flux."  Physically,  this  is  the  convective  he? t flux  which  would  occur 
for  a frozen  boundary  layer  and  a noncatalytic  wall  in  the  absence  of  mass 


transfer;*  it  excludes  all  chemical  energy  contributions.  The  term  q is 

sen 

perhaps  more  usually  written  in  the  form 

*sen  " PeueCH(Hsr  " hsw)  (2_59> 

but,  since  generally  it  is  more  convenient  for  the  user  to  input  Hr  rather  than 

H , q in  Equation  (2-58)  has  been  written  in  a modified  form  in  which  H 
s sen  r 

appears.  This  form  has  the  additional  advantage  that  the  driving  force  for 
energy  transfer  involves  only  edge  gas  states.  The  derivation  of  the  modified 
form  from  Equation  (2-58)  is  given  in  Reference  20  and  Reference  21. 

The  quantity  h^  in  the  qsen  term  is  part  of  the  input  thermochemical 
data  discussed  below.  The  transfer  coefficient  p ti  C_  and  the  recovery  enthalpy 
Hr  are  time  dependent  variables  computed  in  the  program  for  each  analysis  loca- 
tion. The  transfer  coefficient  is  automatically  modified  from  the  nonblown 
value  to  implicitly  account  for  the  effect  of  the  computed  ablation  rates.  The 
following  relation  is  used: 


CH  in(l  + 2aB£c) 
CH,0  2XB^c 


(2-60) 


where 

= implicitly  determined  normalized  thermochemical  ablation  rate 
(“tc/PeueCM) 

= an  input  number  discussed  below 
= ratio  of  blown  to  nonblown  Stanton  number 

Specified  values  of  X allow  the  user  to  fit  blowing  correction  curves 
of  C„/Cu  „ versus  B'  to  account  for  special  effects  in  the  few  cases  where 

n fl|U  cc 

these  are  known  with  confidence,  such  as  molecular  weight  effects  or  variable 
property  effects.  In  view  of  the  uncertainties,  it  is  generally  recommended 
that  X = 0.5  be  used  for  laminar  flow.  A value  X = 0.4  appears  to  correlate 
constant  properties  for  turbulent  data  somewhat  better.  For  graphite  in  air, 
studies  have  indicated  that  a value  of  0.7  for  both  laminar  and  turbulent 
flow  is  most  appropriate. 


More  generally  in  the  presence  of  chemical  reaction  it  is  the  diffusive  heat 
flux  from  the  gas  to  the  wall  even  in  the  presence  of  net  mass  transfer,  pro- 
vided the  boundary  layer  is  frozen  and  the  wall  is  catalytic. 
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The  term  %ond  in  Equation  (2-58)  is  obtained  from  the  in-depth  conduc- 
tion analysis  as  a function  of  Tw  and  S (Equation  (2-54)). 

The  term  qchem  in  Equation  (2-58)  represents  the  net  amount  of  chemical 
energy  fluxes  at  the  surface.  The  z*-dif ference  term  represents  transport  of 
cnemical  energy  associated  with  chemical  reactions  at  the  well  and  ^n  the 
boundary  layer , it  is  the  chemical  energy  parallel  to  the  sensible  convective 
heat  flux  term.  The  z*  driving  forces  for  diffusive  mass  transfer  include  the 
effects  of  unequal  diffusion  coefficients;  for  equal  diffusion  coefficients  the 
reduce  to  the  familiar  mass  fractions  K^.  The  term  represents  energy 

leaving  the  surface  in  the  gross  motion  (blowing)  of  the  gas  adjacent  to  the 
surface.  The  mass  transfer  coefficient  (Pe«eCM)  is  obtained  from  the  blown  heat 
coefficient  ( Peue^jj)  using  a user  specified  factor,  C^/Cy  Remaining  quantities 
not  yet  discussed  are  B^.c,  h^.  r>£ehi  » 2ZziwhiW»  ***  hw*  Tw  does  not  aPPear 
explicitly  but  is  necessary  to  .luate  the  temperature  dependent  values  of 
various  quantities.  The  quantities  Ty,  and  hw  are  input  by  the  user 

as  the  dependent  variables  in  a table  with  three  independent  variables:  p, 

Pe,JeCy  and  B^_  (if  no  chemical  kinetic  effects  are  considered  only  two  inde- 
pendent variables  p and  B^c  are  required). 

T 

Similarly,  the  quantities  and  h^  are  input  as  the  dependent 

variables  in  a table  with  p and  T at  independent  variables.  These  tables  are 
typically  generated  by  the  thermochemistry  codes  described  below.  Further 
discussion  of  these  tables  is  given  in  Section  3.1.9. 

Notice  that  the  erosion  mass  loss  rate  (rae)  does  not  appear  in  Equation 
(2-58) . This  is  because  the  eroded  material  is  assumed  to  leave  the  surface 
with  the  enthalpy  of  the  solid  (hc)  and,  hence,  cancels  with  the  net  incoming 
mass  rate  (m)  to  give 


m h = 
e c 


mtchc 


(2-61i 


where  n>tc  is  only  the  thermochemical  portion  of  the  total  mass  rate  (m) . The 
total  mass  rate  (m)  (and  hence  m^)  is  required,  however,  in  the  conduction  equa- 
tion (Equation  (2-56))  to  compute  the  net  recession  rate  (S) . The  calculation 
of  the  erosion  mass  loss  rate  (me)  is  described  in  Section  2.4.2. 

The  surface  energy  balance  solution  procedure  may  be  summarized  as  fol- 
lows: 

1.  Obtain  Hr,  PeueCH  Q,  and  p from  environment  definition  routines. 

2.  Reduce  the  in-depth  conduction  matrix  to  calculate  the  constants  in 
the  equation  for  qcond  (Equation  (2-54)). 
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3.  Calculate  erosion  mass  loss  rate  (i  ) 

e 

4.  Correct  or  adjust  p u C„  n for  blowing  effects 

6 6 <1|U 

5.  Compute  peueCM  = «VCH) 

6.  Assume  B^c 

7.  With  p,  P^u  C„,  and  B'  , look  up  in  input  surface  thermochemistry 

6 6 N uC 

tables  values  of  Tw#  LziwhiW-  hw 

8.  With  p and  T , look  up  in  input  edge  ga3  thermochemistry  table 
values  of  £ziehiW'  he* 

9.  Construct  Equation  {2-53} , noting  departure  from  zero,  if  any 

10.  Adjust  B^c  guess  to  reduce  departure  from  zero  (Newton-Raphson 
correction) 

11.  Go  to  Step  4 and  continue 

This  procedure  converges  on  a new  B£c  value  in  very  few  iterations.  The  same 
procedure  may  be  used  with  Tw  as  the  independent  variable  and  B^.c  as  a depen- 
dent variable. 


2.4.2  Erosion  Modeling 


Surface  erosion  due  to  hydrometer  particle  impacts  is  modeled  by  two 
different  types  of  correlations  depending  on  the  surface  material.  For  gra- 
phite type  brittle  materials  the  erosion  mass  loss  correlation  is  of  the  form 


A. u^mc 
1 ® p 


sin  0^ 


(2-62) 


where 

m = erosion  mass  loss  flux 
e 

m^n  = p^jsin  0 = incoming  particle  mass  flux 
Pc  = mass  of  particles  per  unit  volume  of  air 
u = vehicle  velocity 

OD 

m = ltd *p  /G  = individual  particle  mass 

p prp 

dp  = particle  diameter 

0 •*  local  body  angle  relative  to  the  axis  (0  = 90°  at  the  stagnation 
point) 
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Tha  constants  (A^ , b,  c and  d)  in  Equation  (2-62)  are  determined  by  correlation 
o£  ground  test  data  for  each  given  material.  For  materials  which  char  (i.e. , 
carbon  phenolic)  a different  set  of  constants  is  required  for  the  charred  and 
virgin  plastic  materials.  The  two  erosion  rates  predicted  for  a fully  charred 
and  virgin  surface  are  "bridged"  together  based  on  the  relative  rates  of  sur- 
face erosion  and  in-depth  char  generation.  Additionally  the  body  angle  (6) 
dependence  in  Equation  (2-62)  does  not  fully  collapse  all  carbon  phenolic  data 
and,  hence,  low  c.nd  high  angle  erosion  correlation  constants  are  required  as 
well  as  a "bridging"  function.  Currently  specific  erosion  correlation  constants 
and  bridging  functions  are  built  into  the  code  for  carbon  phenolic.  Refer- 
ence 19  discusses  the  details  of  this  modeling  and  Reference  23  covers  the 
graphite  erosion  models.  The  input  section  (Section  3.1)  describes  how  these 
specific  correlations  may  be  evoked. 

For  malleable  type  metal  materials  the  erosion  mass  loss  correlation 
is  of  the  form 


where 

Cjj  = deunage  coefficient 

The  damage  coefficient  (Cj,)  is  determined  from  experiment  and  is  primarily  a 
function  of  surface  temperature  (T^) . Typically  the  deunage  coefficient  decreases 
as  the  surface  temperature  approaches  the  melt  temperature.  Built-in  values 
of  CM  versus  Tw  are  available  for  tungsten.  Reference  19  gives  the  details 
of  the  tungsten  correlation  and  the  input  section  tells  how  it  may  be  evoked. 

The  input  required  for  erosion  calculations  is  the  cloud  profile,  which 
is  a table  of: 

• Mass  concentration  (pc) 

• Particle  diameter  (d  ) 

P 

• Particle  specific  gravity  (y#) 
versus  altitude. 

2.4.3  Use  of  Thermochemistry  Codes  to  Generate  Input  Data 

Section  2.4.1  above  makes  it  clear  that  some  complex  tabular  thermo- 
chemical  input  is  required  if  the  surface  energy  balance  boundary  condition  is 
to  be  used.  These  tables  are  generated  by  any  one  of  a number  of  separate 
computer  codes.  The  most  recent  such  code  is  designated  General  Nonequilibrium 
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Ablation  Thermochemistry  Code  (GNAT).  Tt  is  a general  open  and  closed  system 
thermochemical  nonequilibrium  code  specifically  constructed  for  this  purpose 
(Reference  24).  Other  Aerothcrm  thermochemistry  codes  which  treat  only  equi- 
librium thermochemistry  are  in  existence.  The  most  recent  o.  these  is  the 
Equilibrium  Surface  Thermochemistry  Code,  Version  3 (EST3) , which  is  described 
in  Reference  25.  A generally  similar  code  which  differs  from  EST3  only  in  added 
detail  is  designated  ACE  and  is  described  in  Reference  26.  An  older  version  of 
EST3  was  designated  EST2  and  is  described  in  Reference  27.  To  obtain  the  nec- 
essary data  tables  for  input,  the  user  selects  sets  of  values  for  the  pressure 
(P) , transfer  coefficient  (PeueCM),  and  nondimensional  thermochemical  ablation 
rate  (B^c) . (Note,  if  chemical  kinetics  are  not  considered  the  only  parameters 
required  are  pressure  (p)  and  ablation  rate  (B^J.) 

The  user  specifies  the  elemental  composition  of  the  environment  gas  and 
the  ablating  material,  and  supplies  some  general  species  thermochemical  data 
for  all  molecules  to  be  considered  in  the  system.  Finally,  the  user  specifies 
the  unequal  diffusion  coefficients  if  they  are  important.  The  thermochemistry 
code  then  computes  all  the  dependent  quantities  of  interest  at  each  table  point 
in  the  p x B ' matrix  of  independent  variable  values,  namely,  T , h.w,  and 

tC  W ^ 

hw,  and  punches  this  information  on  cards.  Similarly,  the  tables  of  2lziehiW 
and  h^  values  are  prepared  as  functions  of  p and  T,  and  punched  out  on  cards. 

All  these  cards  form  part  of  the  Table  09  card  input  deck  (see  Section  3.1.9). 

2.4.4  Simpler  Forms  of  the  Surface  Energy  Balance  Equation 

As  noted  in  Section  2.5.2  above,  for  equal  diffusion  coefficients  the 

zf  driving  forces  reduce  to  the  simple  mass  fractions  K. . If  in  addition  to 
X x i* 

equal  diffusion  the  user  specifies  that  p u C„  = p u c„,  then  since  T'K.  h.w  = h» 

«p  c € n c 6 M 1C  X w 

and  £Kiwb^s  = hw  by  definition.  Equation  (2-58)  simplifies  to  the  more  familiar 
form 

PeueCH  (Hr  ‘ (1  + Btc}hw)  + \chc  “ <*cond  + °wqrad  " FoeTw  = 0 (2‘64) 

T 

In  this  expression  h^  and  Sz*eh£W  do  not  appear,  hence  the  corresponding 
table  is  not  necessary  and  need  not  be  included  in  the  input  (see  Section  3.1.9 
below) . 

A steady  state  ablation  option  is  also  available.  If  this  option  is 
specified  the  <3concj  term  in  Equation  (2-58)  is  calculated  by  taking  an  energy 
balance  on  a control  volume  extending  from  below  the  ablating  surface  down  to 
the  thermally  unaffected  material  (see  the  following  sketch) . 
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Sketch  of  Control  Volume  Around  Thermally  Effected  Material 


Energy  conservation  on  the  above  control  volume  gives 


*cond 


dE 

dt 


(2-65) 


where 

Ti 

h^  » the  enthalpy  of  the  thermally  uneffected  material  before 
exposure  (T^  is  assumed  to  be  530°R) 

~ rate  of  energy  storage  in  control  volume 

The  steady  state  assumption  implies  that  dE/dt  is  zero  and  corresponds  to  the 
physical  situation  when  the  temperature  profile  relative  to  the  moving  surface 
is  invariant  with  time.  The  assumption  is  accurate  for  low  conductivity  abla- 
tors and  for  high  ablation  rate  situations.  By  considering  dE/dt  * 0,  <lcon(j 
may^be  calculated  from  Equation  (2-b5)  as 


‘cond 


“ (hc  - hc0 


(2-66) 


Notice  that  for  the  steady  state  assumption,  <3con<j  is  independent  of  material 
thermal  properties  and  response  history. 

2.5  SURFACE  POINT  MOVEMENT  AND  SURFACE  SMOOTH Ihd 

The  surface  energy  balance  determines  the  recession  rate  normal  to  ^he 
surface  at  each  body  calculate-.,  point.  Based  on  the  time  step  size  (At)  these 
points  are  then  moved  the  _ espond  .ng  distance  to  define  new  body  points  at 
the  end  of  the  time  step  ” r-  : c-*-  jcdy  points  are  then  used  to  define  new  sur- 
face inclination  angle'  he  body  points.  In  a typical  nosetip  shape  change 
problem  the  size  of  important  geometric  features  in  the  stagnation  region  de- 
creases in  turbulent  flow,  and  eventually  becomes  smaller  than  can  be  efficiently 
modeled  with  typical  body  point  spacing.  When  a numerically  limited  nose  radius 
is  reached  logic  is  applied  to  define  an  apparent  nose  radius. 


The  numerics  of  shape  change  are  described  in  the  following  sections. 

The  shape  change  geometry  for  body  point  movement  is  described  in  Section  2.5.1; 
surface  angle  definition  techniques  are  presented  in  Section  2.5.2;  and  the 
apparent  nose  logic  is  discussed  in  Section  2.5.3. 

2.5.1  Shape  Change  Geometry 

With  reference  to  Figure  2-6,  the  location  and  shape  of  the  surface  is 
completely  defined  by  A(X,$,t).  The  rate  of  change  of  A with  time  can  be  re- 
lated to  the  surface  normal  recession  rate,  S by  the  following  equation: 


3A 

at 


aA/axV  ( 3A/a»\2 
( l+AA  j yR+AcosB/ 


(2-67) 


where  A,  R and  6 are  the  geometric  parameters  of  the  internal  contotr  defined, 
respectively,  as:  curvature,  radial  distance  from  the  nosetip  axis  end  angle 

of  inclination  with  respect  to  the  nosetip  axis. 

The  Equation  (2-67)  is  written  in  an  explicit  finite-difference  form 
and  solved  for  the  values  of  A at  the  n+1  time  step.  Along  the  centerline 
where  R » 0,  Equation  (2-67)  is  not  applied  and  the  condition  3A/3X  = 0 is 
used  to  determine  the  new  position  of  the  stagnation  point. 

In  the  steady  state  conduction  option  the  body  points  are  moved  along 
lines  of  constant  radius  as  indicated  by  the  sketch  below. 

e 


<L 


- z 


The  relation  for  the  amount  of  axial  movement  is 


normal 
sin  6 


As  “ 


At 


(2-68) 


When  using  the  steady  state  energy  balance,  conduction  considerations  do 
not  limit  time  step  size;  hence,  the  only  consideration  limiting  time  step  size 
is  shape  change.  In  other  words,  the  calculated  recession  rate  distribution 
cannot  be  applied  over  such  a time  span  that  the  body  shape  (and,  hence,  reces- 
sion rate  distribution)  changes  in  a drastic  manner.  The  criteria  applied  is 
that  the  tangent  of  the  local  body  angle  may  not  change  by  more  than  a factor 
of  two. 


2.5.2  Surface  Angle  Definition 

Numerical  shape  change  calculation  techniques  are  strongly  sensitive  to 
the  method  used  to  define  the  local  surface  angle  since  this  angle  strongly 
influences  the  surface  pressure,  heat  transfer,  ablation,  and  erosion  calcula- 
tions. Circular  curve  fits  ar.d  straight  line  interpolation  techniques  are 
available. 

The  circular  curve-fit  method  involves  basically  fitting  a circular  arc 
through  the  point  of  interest  and  the  points  on  either  side  (i.e.,  three  points 
define  a circle) . The  body  angle  is  then  defined  as  the  tangent  to  the  circle 
at  that  point.  Exceptions  are  taken  to  this  definition  if  the  radius  of  curva- 
ture is  negative  in  order  to  avoid  unrealistic  concave  shapes. 

2.5.3  Apparent  Nose  Model 

As  shape  change  proceeds  on  nosetips  for  which  transition  is  near  the 
nose,  the  stagnation  point  radius  of  curvature  becomes  too  small  to  model  with 
the  typical  body  point  spacing.  The  code  has  internal  logic  to  determine  when 
this  numerically  limited  nose  radius  is  reached,  and  at  that  time  an  effective 
spherical  nose  radius  is  computed.  This  specification  controls  only  the  detail 
at  the  stagnation  point  and  does  not  limit  or  redefine  the  overall  shape  of 
the  nosetip. 

The  apparent  nose  radius  logic  used  with  the  circular  curve  fits  involves 
basically  fitting  a tangent  sphere  into  the  "cone"  formed  by  the  second  and 
third  body  points,  and  is  primarily  based  on  geometrical  considerations. 

Neither  of  the  apparent  nose  radius  or  body  angle  definition  techniques 
described  above  is  completely  successful  in  predicting  all  observed  nosetip 
shape  change  regimes.  Hence,  both  must  still  be  considered  to  be  in  the  de- 
velopmental stage. 
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SECTION  3 

DESCRIPTIONS  OF  INPOT  AND  OUTPUT 

This  section  provides  detailed  user  oriented  instructions  for  code  input 
and  a description  of  the  output.  The  input  instructions  are  described  in  Sec- 
tion 3.1  and  output  features  are  covered  in  Section  3.2. 

3.1  INPUT  INSTRUCTIONS 

The  input  to  the  code  can  be  read  either  from  data  cards  for  an  initial 
run  or  from  magnetic  tape  or  disk  for  a restart  run.  The  details  of  the  input 
for  each  of  these  types  of  runs  are  described  below.  The  basic  input  for  an 
initial  run  ronsists  of: 

• One  restart  information  card. 

• Three  title  cards. 

• Nine  input  tables. 

Not  all  nine  of  the  input  tables  are  required  for  every  run.  Each  table  is 
preceded  by  a single  card  containing  the  identifying  table  number. 

For  a restart  run  only  the  single  restart  information  card  is  required, 
the  rest  of  the  information  is  read  from  magnetic  tape  or  disk. 

The  following  sections  describe  the  restart  information,  title  cards, 
and  nine  input  tables,  respectively. 

3.1.1  Restart  Information 

The  "restart'"  card  is  the  first  card  in  the  data  deck.  If  the  run  is  a 
restart  it  is  the  only  card  required,  and  tells  the  code  the  iteration  from 
which  to  begin  restart.  For  an  initial  run  the  restart  card  tells  the  code  how 
often  to  write  restart  files.  All  restart  reading  and/or  writing  is  done  on 
Logical  Unit  11  and  it  should  be  assigned  accordingly.  The  restart  card  for- 
mat, is  as  follows: 
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Card 


NO. 

Columns 

Format 

Data 

Units 

1 

1-3 

13 

ISKI’  - number  of  environmental  calls 
between  restart  file  writes. 

— 

4-6 

13 

IRESRT  - restart  flag 

— 

0 - first  run  no  restart 

>0  - transient  restart  - read  the 

IRESRTth  set  of  data  on  Unit  11 
to  start  the  run.  No  other  in- 
put is  required 

<0  - steady  state  restart 

Note  that  if  both  ISKIP  and  XRESRT  are  equal  to  zero  no  data  will  be  read  or 
written  on  Unit  11  and  it  need  not  be  assigned. 

3.1.2  Title  and  Heading  Information 

The  second  set  of  input  data  are  three  title  cards.  They  are  used  to 
transmit  title  and  heading  information  to  the  output.  The  first  72  columns  of 
each  of  these  cards  may  be  used  for  the  title,  the  alphameric  information  in 
columns  61  through  72  of  the  third  card  being  used  as  a page  heading  on  all 
pages  after  the  first. 

3.1.3  Table  01  - General  Program  Constants 

These  cards  supply  the  code  with  computation  time  information  and  program 
flags  which  indicate  options  to  be  subsequently  read. 


Card 


Ko. 

Columns 

Format 

Data 

Units 

1 

1-2 

12 

Table  No. 

— 

2 

1-12 

E12.5 

Initial  value  of  problem  time 

sec 

13-24 

E12.5 

Final  value  of  problem  time 

sec 

— 

25-36 

E12.5 

First  output  tine  increment.  This  in- 
terval represents  the  time  increment 
for  output.  Provision  for  changing 
this  time  increment  within  a run  is 
provided  by  the  NTIC  flag  described 
below. 

37-48 

E12.5 

DLTMIN  - rime  step  stability  flag 

sec 

<0  - no  stability 

=0  - set  to  10" * sec 

>0  - stability  criteria  used 
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Table  01  (continued) 

Card 

Mo.  Columns  Format 
49-60  E12.5 

61-72  E12.5 


3 1-3  13 


4-6  13 


Data 


Units 


CTF  - defined  by  Equation  (2-56). 

If  CTF  is  less  than  1.2  or  greater 
them  1.7  it  is  set  to  1.3. 


STRD  - maximum  desired  surface  tem- 
perature rise  in  one  conduction  time  °R 
step.  Zf  STRD  is  less  than  49  or 
greater  than  201  it  is  set  to  75°R. 

STRD  is  not  required  for  the  steady 
state  option. 


TC  - Flag  denoting  type  of  transi- 
tion criterion  to  be  subsequently 
input  in  Table  04.  TC  < 0 denotes 
transitional  heating  is  used  in  the 
surface  energy  balance  and  TC  > 0 
denotes  abrupt  transition  is  used. 

ABS  (TC) 


0 

.1 

2 

3 

4 


5 


6 


all  laminar  flow  (Table  04  is 
not  needed) 

momentum  thickness  Reynolds 
No.  vs.  edge  Mach  No. 

run  length  Reynolds  No.  vs. 
edge  Mach  No. 

axial  distance  vs.  altitude 

rough  wall  transition  based 
on 

/ s \ ^ ( 2300,  onset 

Rek  \5*/  * 12000,  location 

(Table  04  is  not  needed) 

rough  wall  transition  based  on 

0.7 

m(255,  onset 
(215,  location 

(Table  04  is  not  needed) 

fully  turbulent  flow  (Table  04 
is  not  needed 


Re£ 


(* i + i + §1  -s)  0 

_\10  « dJ  _ 


ENV  - flag  denoting  environment  option 
to  be  subsequently  read  in  Table  02 

1 - flight  option 

2 - wind  tunnel 
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Table  01  (continued) 

Card 

No.  Columns  Format 


7-9  13 


10-12  13 


13-15  13 


16-18  13 


19-21  13 


Data 


Units 


3 - ballistic  range 

4 - general 

5 - arc  heater 

CF  - flag  controlling  curve  fit  and 
apparent  nose  option 

0 - circular  curve  fits  and  no 
apparent  nose  logic 

2 - circular  curve  fits  and  appar- 
ent nose 

SO  - special  output  flag 

0 - boundary  layer  solution  only 

(no  ablation  or  shape  change) . 

1 - general  problem  with  ablation 

and  shape  change.  Shape  pro- 
files written  on  file  15 

NT1C  - number  of  time  interval  changes 
"(not  number  of  time  intervals) 

NTIC^jj  * 10.  A non-zero  entry  in  this 
column  causes  sets  of  time  interval 
changes  to  be  read  from  the  next  card. 

1SS  - conduction  option  flag 

0 - transient  conduction  option. 

Sphere-cone  initial  geometry 
with  geometric  orogression  dis- 
tributions of  surface  and 
in-depth  grids. 

1 - steady  state  conduction  option. 

Initial  geometry  and  surface 
point  distributions  same  as  above. 

2 - steady  state  conduction  option. 

General  initial  geometry  and 
surface  points  distribution.  The 
details  of  this  input  are  de- 
scribed in  Section  3.1.5. 

IPRNT  - flag  which  determines  the  amount 
of  environmental  output  at  print  times. 
Six  output  tables  are  available  and  the 
contents  of  each  is  described  in  Sec- 
tion 3. 2. 2.1.  IPRNT  < 0 denotes  output 
for  each  integration  point  and  IPRNT  > 0 
denotes  output  for  body  points  only. 


3-4 


4 - print  Tables  1,  3,  4,  5 

and  6. 


22-24  13  LPRNT  - flag  similar  to  IPKNT  which  de-  — 1 

termines  the  amount  of  output  at  inter-  j 

mediate  computation  times  (i.e. , when  ,1 

' confutation  time  not  equal  to  print  time) . 3 

LPRNT  < 0 denotes  output  for  each  inte-  | 

gration  point  and  LPRNT  > 0 denotes  out- 
put for  body  points  only. 

ABS  (LPRNT) 

0 - no  output 

1 - abbreviated  output  of  environ- 

ment and  recession  only 

2 - print  Tables  1 and  2 

3 - print  Tables  1,  3,  4 and  5 

4 - print  Tables  1,  3,  4,  5 

and  6. 

4 This  card  supplies  information  for  changes  in  the  output  time  inter- 

val and  is  read  only  if  NTIC  >0. 


1-12 

E12.8 

Second  output  time  interval 

sec 

13-24 

E12.8 

Time  for  change  to  second  output  time 
interval 

sec 

25-36 

E12.8 

Third  output  time  interval 

sec 

37-48 

E12.8 

Time  for  change  to  third  output  time 
interval 

sec 

49-60 

E12.8 

Fourth  output  time  interval 

sec 

61-72 

E12.8 

Time  for  change  to  fourth  output  time 
interval 

sec 

5 Same  for  NTIC  output  time  interval  changes.  Three  sets  per  card,  to 

(etc)  a maximum  of  ten  sets. 


3.1.4  Table  02  - Environment  Table 

The  basic  environment  information  required  by  the  code  is  the  freestream 
state  (pressure  and  density'  and  vehicle/gas  relative  velocity.  Given  this  \ 
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Table  02  (continued) 

information  the  code  performs  real  gas  calculations  for  air  to  find  the  stag- 
nation conditions.  To  aid  the  user  in  performing  calculations  for  various 
common  flight  and  ground  test  facility  environments,  five  environment  input 
options  are  provided.  They  are: 

1.  Flight  environment  - input  altitude  and  velocity  as  a function  of 
time;  free  stream  conditions  a; e found  from  a built-in  ARDC  standard 
atmosphere  table. 

2.  Hind  tunnel  environment  - input  supply  pressure  and  temperature  as  a 
function  of  time  as  well  as  free  stream  Mach  No.  and  ratio  of  specific 
heats  (y) ; free  stream  conditions  are  found  based  on  assumed  air 
isentropic  expansion  at  constant  y. 

3.  Ballistic  range  environment  - input  projectile  velocity  and  range 
pressure  as  a function  of  time;  air  assumed  to  be  at  75°F  to  obtain 
free  stream  density. 

4.  General  environment  - input  free  stream  pressure , density,  and  veloc- 
ity as  a function  of  time. 

5.  Arc  heater  environment  - input  consists  of  a quantitative  description 
of  the  arc  heater  flow  field  and  of  the  start-up  transient  which  the 
model  experiences  when  injected  into  this  flow  field. 

Tor  the  flight,  wind  tunnel,  ballistic  range,  and  general  environment 
options.  Table  02  is  basically  a time  table  of  environment  conditions.  For 
the  arc  heater  option  the  environment  is  more  complex  and  in  this  case  Table  02 
provides  input  of  both  fch2  spacial  and  temporal  variation  of  environment. 

There  are  a maximum  of  50  entries  in  this  table;  at  least  two  entries 
are  required. 

3. 1.4.1  Input  for  Flight  Option,  ENV  * 1 
Card 


No. 

Columns 

Format 

pata 

Units 

1 

1-2 

X2 

Table  No.  02 

— 

2 

1-2 

12 

Must  be  blank 

— 

3-14 

E12.8 

Time 

sec 

15-26 

E12.8 

Altitude 

ft 

27-38 

E12.8 

Velocity 

ft/sec 

3 

Same  as 

Card  No.  2 

for  increasing  time. 

— 

(etc) 
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Table  02  (continued) 

3. 1.4. 2 Input  for  Wind  Tunnel  Option,  ENV  - 2 
Card 


No. 

Columns 

Format 

Data 

Units 

1 

1-2 

12 

Table  No.  02 

— 

2 

1-2 

12 

Must  be  blank 

— 

3-14 

E12.8 

* Time 

sec 

15-26 

E12.8 

Supply  pressure 

psia 

27-38 

E12.8 

Supply  temperature 

•F 

39-50 

E12.8 

Free  stream  ratio  of  specific  heats 
(Card  No.  2 only) 

— 

51-62 

B12.8 

Free  stream  Mach  No.  (Card  No.  2 only) 

— 

3 

1-2 

12 

Must  be  blank 

— 

3-14 

E12.8 

Tine 

sec 

15-26 

E12.8 

Supply  pressure 

psia 

27-28 

E12.8 

Supply  temperature 

"F 

4 

tetc) 

Same  as 

Card  No . 3 

for  increasing  time. 

3. 1.4. 3 

Input  for  Ballistic  Range  Option,  ENV  = 3 

Card 

No. 

Columns 

Format 

Data 

Units 

1 

1-2 

12 

Table  No.  02 

— 

2 

1-2 

12 

Must  be  blank 

— 

3-14 

E12.8 

Time 

sec 

15^26 

E12.  8 

Range  static  pressure 

atm 

27-38 

E12.8 

Projectile  velocity 

ft/sec 

3 Sane  as  Card  No.  2 for  increasing  tine. 

(etc) 

3. 1.4. 4 Input  for  General  Environment  Option,  ENV  = 4 
Card 


No. 

Columns 

Format 

Data 

Units 

1 

1-2 

12 

Table  No.  02 

— 

1-2 

12 

Must  be  blank 

— 

3-7 
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Table  02 

(continued) 

Card 

No. 

Columns 

Format 

Data 

Units 

3-14 

E12.8 

Time 

sec 

15  26 

E12.8 

Free  stream  static  pressure 

atm 

27-38 

E12.8 

Free  stream  static  density 

lb/ft* 

39-50 

E12.8 

Free  stream  velocity 

ft/sec 

51-62 

E12.8 

Ratio  of  specific  heats  (Card  No.  2 
only) 

— 

3 

1-2 

12 

Must  be  blank 

— 

3-14 

E12.8 

Time 

sec 

15-26 

E12.8 

Free  stream  static  pressure 

atm 

27-38 

E12.8 

Free  stream  static  density 

lb/ft  * 

39-50 

E12.8 

Free  stream  velocity 

ft/sec 

4 

(etc) 

Same  as  Card  No.  3 

for  increasing  time,  maximum  of  50  entries. 

3. 1.4. 5 

Input  for  Arc  Heater  Environment 

The  free  stream  environment  produced  by  an  arc  plasma  generator,  such  as 
the  AFFDL  50  MW  RENT  facility,  has  characteristics  distinctly  different  from  the 
environment  options  described  above  (ENV  = 1-4) . In  general,  the  pressure  level 
is  constant  with  time  but  varies  with  distance  from  the  nozzle  exit  because  of 
nonparallel  flow  streamlines.  Input  consists  of  one  card  containing  the  steady 
operating  conditions  (including  model  location)  and  a table  of  normal  shock  total 
pressure  ratio  as  a function  of  distance  from  the  nozzle  exit. 

An  option  also  exists  for  specifying  a free  stream  pressure  variation 
during  a start-up  transient.  The  ratio  of  instantaneous  total  free  stream  pres- 
sure to  the  steady ‘value  for  two  or  more  times  are  input.  The  option  is  flagged 
by  reading  in  a nonzero  entry  for  the  length  of  the  start-up  transient  (DTIME) 
which  is  read  from  the  second  card.  The  values  of  total  pressure  ratio  versus 
time  are  input  using  the  same  read  statement  which  specifies  the  spacial  varia- 
tion of  pressure  ratio,  by  adding  100  to  each  of  the  times.  Again  a maximum  of 
50  entries  is  allowed. 
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Table  02  (concluded) 
Card 


No. 

Columns 

Format 

Data 

Units 

1 

1-2 

12 

Table  No.  02 

— 

2 

1-2 

12 

Must  be  blank 

— 

3-14 

E12.8 

Pt«ateady'  tota*-  *ree  stream  pressure 

atm 

15-26 

E12.8 

Hq,  total  enthalpy 

Btu/lb 

27-38 

E12.8 

Xq,  axial  distance  from  nozzle  exit 
plane  to  stagnation  point  (assumes 
model  advance  with  recession) 

in. 

39-50 

E12.8 

Ki,  multiplying  constant  in  the  Fay  & 
Riddell  stagnation  point  heating  rela- 
tion (if  Ki  * 0,  it  is  set  to  1.0) 

51-62 

E12.8 

DTIME  - Duration  of  start-up  transient 

sec 

63-74 

E12.8 

Y,  specific  heat  ratio  (if  y = 0,  it 
is  set  to  1.2) 

— 

3 

1-2 

12 

Must  be  blank 

— 

3-14 

E12.8 

X,  distance  from  nozzle  exit  plane 

in. 

15-26 

E12.8 

pt2/pt-steady'  pressure  ratio 

— 

4 

(etc) 

Same  as 

Card  No.  3 

for  increasing  X values. 

n 

1-2 

12 

Must  be  blank 

— 

3-14 

E12.8 

Time  + 100 

sec 

. 

15-26 

E12.8 

pt»/pt«8teady/  pressure  ratio 

— 

n+1  Same  as  Card  No.  n for  increasing  time,  maximum  of  50  entries  includ- 
ing Card  No.  3.  There  must  be  at  least  two  x-values  and  two  time 
values  when  this  option  is  used.  Time  variant  environment  logic  is 
not  debugged. 


3.1.5  Table  03  - Geometrv 


In  the  first  part  of  this  section  initial  geometry  and  in-depth  grid  arc 
described  and  in  the  second  part  the  input  format  of  Table  03  is  explained. 


3. 1.5.1  Initial  Geometry 

The  initial  body  geometry  is  input  as  a table  of  coordinates  for  several 
points  on  the  body.  This  input  is  different  for  steady  state  and  transient 
conduction  options.  The  differences  are  described  below. 
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Table  03  (continued) 

1.  Transient  conduction  option  - a sphere  cone  initial  geometry  and 
geometric  progression  distribution  of  body  points  are  assumed. 

Input  sphere  radius,  cone  half  angle,  number  of  body  points  and 
the  common  ratio  of  the  geometric  progression.  The  initial  geom- 
etry may  be  of  a shell  or  plug  category.  For  a shell,  the  internal 
contour  is  assumed  to  be  a sphere  cone  whose  sphere  radius  and  cone 
half  angle  are  input.  For  a plug,  the  sha<r>k  geometry  is  xriput. 

The  conduction  option  flag  in  this  case  is  ISS  * 0.  The  details  of 
these  inputs  are  described  in  Section  3. 1.5.4. 

2.  Steady  state  option  - under  this  option,  three  types  of  input  are 
possible  for  initial  geometry  and  surface  point  distributions: 

• Sphere  cone  geometry  and  geometric  progression  distribution  of 
surface  points.  This  input  is  the  same  as  the  transient  option 
input.  The  conduction  option  flag  in  this  case  is  ISS  = 1. 

• Sphere  cone  geometry  with  uniform  surface  point  distribution. 
Input  sphere  radius,  cone  half  angle,  axial  length  and  number 
of  body  points  desired  on  sphere  and  cone. 

• General  body  - input  a table  of  up  to  30  body  coordinates  (r,z). 
The  conduction  option  flag  in  the  above  two  cases  is  ISS  <=  2. 

3 . 1 . 5 . 2 In-Depth  Grid 

The  input  of  the  in-depth  grid  applies  only  to  the  transient  conduction 
option  of  the  code.  The  in-depth  grid  system  is  shown  on  Figure  2-6.  Since 
the  code  is  presently  limited  to  axisymmetric  geometries,  we  will  only  consider 
the  grid  distributions  in  the  X-  and  n-directions.*  Geometric  progression  dis- 
tributions of  grids  in  both  X-  and  n-directions  are  assumed.  Therefore  it  is 
only  required  to  input  the  number  of  the  grid  points  and  the  common  ratios  in 
X-  and  n**directic.ns. 

3. 1.5. 3 Surface  Temperature 

The  code  assigns  one  input  value  of  surface  temperature  to  all  body 
points  for  all  input  options  except  the  general  body.  If  no  surface  tem- 
perature is  specified,  a default  of  530°R  is  utilized.  For  the  general 
shape  option  a surface  temperature  distribution  can  be  input. 

% defined  on  Page  2-28. 
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Table  03  (continued) 

3. 1.5. 4 Table  03  Input  Format 
Card 

No.  Columns  Format  Data 

1 1-2  12  Table  No.  03 

Read  The  Following  Cards  Only  if  ISS  ° 0 or  1.  If  ISS  = 2 go  to  page  3 


2 1-4  4A  PS  - = "PLUG"  or  "SHEL", 

specifying  the  nosetip  geometry 

9-13  4A  SB  - ignored  if  PS  = "PLUG",  must 

Be  "SOLI"  if  a solid  body  is 
desired.  The  solid  body  is  a shell 
with  no  internal  contour. 

3 1-5  15  ABLATE  - integer  variable  defining 

surface  movement. 


“ 0 surface  movement  is  specified 
in  subroutine  TRANS 

* 1 surface  movement  is  calculated 
in  subroutine  ABL8 

6-10  15  MOVE  - integer  variable  defining 

surface  movement 


«=  0 no  surface  movement 

■ 1 allow  surface  movement 

If  ABLATE  = 1,  MOVE  must  = 1 

11-15  15  KAPFLG  - has  meaning  for  shell 

geometry  only 

= 1 zero  curvature  of  the  inter- 
nal contour  and  BETA  = ir/2, 

ZJ  = XJ,  ZB  = 0 and  KAPPA  = 0 

• 0 finite  curvature  of  the 
internal  contour 


16-20  15  IPHI  - flag  for  extrapolation  in 

^-planes.  Applies  only  to  three- 
• dimensional  in-depth  computations. 

= 0 all  back  extrapolations  to 

be  done  for  half  the  ^-planes; 
the  other  half  are  defined 
symmetrically 

*=  1 the  back  extrapolation  to  be 
done  for  L=1  only;  the  other 
planes  are  set  to  the  L=1 
values 


Units 


14. 
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Table  03  (continued) 

Card 

No.  Columns  Format 


Data 


Units 


4 


5 


21-25 


26-30 


1-5 


15 


15 


15 


NOSYM  - flag  for  symmetry  in  ^ -direction 
Applies  only  to  three-dimensional  shape 
change  calculations. 

= 0 assumes  a symmetry  plane 

■ 1 implies  no  symmetry 

I ZERO  - flag  for  angle  of  attack 

= 1 zero  angle  of  attack; 

T7FLG  and  T22FLG  will  be 
turned  on,  NOSYM  set  to 
zero  and  LMAX-3. 

In  the  present  axisymmetric  calculations, 
a value  of  1 should  bt  entered  for  1ZERO. 

JMAX  - number  of  X steps  +1  (maximum 

sr?3) 


6-10 


15  KMAX  - number  of  n steps  +1  (maximum 

oF76) 


11-15 

15 

LMAX  - number 
oFT) 

of  i 

> steps  +1  (maximum 

1-10 

F10.0 

AE  - controls 
n-direction 

the 

mesh  spacing  in  the 

— 

11-20 

F10.0 

AX  - controls 
SPdirection 

the 

mesh  spacing  in  the 

— 

AE  and  AX  are 

the 

common  ratios  of  the 

geometric  progressions  for  grid  dis- 
tributions in  n-  and  X-directions. 

For  uniform  raosh,  set  AE=AX=1.  For  a 
fine  body  point  distribution  near  the 
stagnation  point  enter  AX>1,  and  for  a 
fine  in-depth  grid  near  the  nosetip 
surface  enter  AE<1. 


Geometric  Parameters  Relative  To  The  Internal  Contour 
(For  a plug  these  variables  are  calculated  internahly 
and  the  input  values  ignored.) 


1-10 

F10.0 

RN  - nose  radius 

ft 

11-20 

F10.0 

THC  - cone  half  angle 

deg 

21-30 

F10.0 

XLEN  - maximum  X-distance.  In  the 

ft. 

&6LID  case,  XLEN  must  be  an  angle 

rad 

in  radians. 
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Table  03  (continued) 

Card 

No.  Columns  Format 

31-40  F10.0 

7 1-10  F10.0 


11-20  F10.0 


Data 


BSUBB  - square  of  the  ratio  of  the 
major  to  minor  axis  of  an  ellipse. 
Enter  a value  of  1 for  spherical 
contour. 

T2FLG  - flag  for  X-variation 

=1  no  X-variation,  TERM2  and 
T2  are  internally  set  to  0 

«*  0 allow  X-variation 

T7FLG  - flag  for  $ -variation 

*=  1 no  ♦-variation,  TERM7  and 
T2  are  internally  set  to  0 

* 0 allow  ♦-variation 

In  the  present  code  enter  a value 
of  1 for  this  flag. 


Cards  8,  9 

i and  10  Are 

To  Be 

Input  For  The  Plug  Geometry  Only 

8 

1-10 

F10. 0 

GAMMA  - angle  from  the  horizontal 
that  defines  the  inclination  of  the 
plug  shank 

11-20 

F10.0 

RLEN  - radius  of  the  plug  shank 

(rx) 

21-30 

F10.0 

ZLEN  - length  of  the  p’ug  shank 

(zx) 

9 

1-5 

15 

KPMAX  - number  of  Z steps  +1 

10 

1-10 

F10.0 

T22FLG  - flag  for  ♦-variation  in 
the  shank 

=1  no  ♦-variation;  T22  is 

equal  to  zero 

set 

= 0 allow  ♦-variation 

11 

1-10 

F10.0 

TINITL  - initial  temperature  of 
body 

the 

Geometric 

Parameters 

Relative  To  The  External  Contour 

12 

1-10 

F10.0 

RN2  - nose  radius 

11-20 

F10.0 

THETA2  - cone  half  angle 

21-30 

F10.0 

2L  - nosetip  overhang  for  SHEL; 
axial  length  for  PLUG  or  SOLID  ( 
Figure  2-6) 

total 

[see 

13 

1-5 

15 

NNMAT  - material  index  assigned 

to 

all  surface  and  in-depth  grid  points 


Units 


deg 

ft 

ft 


°R 

ft 

deg 

ft 
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Table  03  (concluded) 
Card 


No. 

Columns 

Format 

Data 

Units 

Read 

The  Following 

Cards  Only  if  1SS  = 2 

2 

1-5 

15 

NS  - number  of  points  on  the  body 
surface  (maximum  30  points) 

— 

>0  sphere-cc.ne  geometry 

40  general  shape 

6-10 

15 

NPN  - number  of  points  on  the  nose; 
applicable  only  to  sphere-cone  option 
(NS>0) 

11-20 

F10. 5 

RSTAGZ  - initial  nose  radius 

ft 

21-30 

F10.5 

ZMAX  - maximum  axial  length 
(sphere -cone  option  only) 

ft 

31-40 

F10. 5 

ANGLI  - initial  cone  half  angle 
(sphere-cone  option  only) 

deg 

41-50 

F10.5 

TS  - initial  body  temperatures 
l7  entered  zero,  it  is  set  to  530°R 

°R 

General  Shape  Option  Only  - 

(read  only  if  NS40) 

3 

1-2 

12 

NC  - flag  to  read  the  coordinates  of 
tKe  body  points 

= 0 keep  reading 

/ 0 stop  reading.  This  indicates 

that  the  card  is  the  last  of 
its  kind. 

3-14 

E12.8 

ZSP  - body  point  axial  length  measured 
from  the  stagnation  point 

ft 

15-26 

E12.8 

RSP  - body  point  radial  length 

ft 

27-38 

E12.8 

ATS  - body  point  temperature.  If 
entered  zero,  it  will  be  set  to  TS 

°R 

39-40 

12 

I MAT  - body  point  material  index 

— 

4 

Same  as  Card  No.  3 

for  the  rest  of  the  body  points. 

(etc) 


3.1.6  Table  04  - Transition  Criteria 


This  table  communicates  criterion  for  specifying  when  boundary  layer 
transition  occurs.  Of  the  six  transition  criteria  available  three  require  tab- 
ular information  which  is  transmitted  through  this  table;  they  are  momentum 
thickness  Reynolds  No.  (Reg)  vs  local  Mach  No.  (Mg) , run  length  Reynolds  No. 
(Re_)  vs  local  Mach  No.  (M  ) and  axial  transition  location  vs  altitude.  The 

S 6 

TC  flag  read  with  the  general  constants  denotes  which  criterion  applies.  This 
table  must  contain  at  least  two,  but  no  more  than  30  entries. 


3. 1.6.1 

Input  for 

Reg  versus 

Me  (ABS (TC)  = 1) 

Card 

No. 

Columns 

Format 

Data 

Units 

1 

1-2 

12 

Table  No.  04 

— 

2 

1-2 

12 

Must  be  blank 

— 

3-14 

E12.8 

Local  edge  Mach  No.  (Me) 

— 

15-26 

E12.8 

Transitional  momentum  thickness 
Reynolds  No.  (Re0) 

— 

3 

(etc) 

Same  as  Card  No.  2 for  increasing  Mg. 

3. 1.6.2 

Input  for  Rea  versus 

Me  (ABS (TC)  - 2) 

Card 

No. 

Columns 

Format 

Data 

Units 

1 

1-2 

12 

Table  No.  04 

— 

2 

1-2 

12 

Must  be  blank 

— 

3-14 

E12.8 

Local  edge  Mach  No.  (Me> 

— 

15-16 

E12.8 

Transitional  run- length  Reynolds  No. 

— 

(Res) 

3 Same  as  Card  No.  2 for  increasing  M . 

e 


3. 1.6. 3 Input  for  Axial  Transition  Location  vs  Altitude  (ABS(TC)  * 3) 
Card 


No. 

Columns 

Format 

Data 

Units 

1 

1-2 

12 

Table  No.  04 

2 


1-2 


12 


Must  be  blank 
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Table  04  (concluded) 

Card 

No.  Columns  Format 


Data 


Units 


3 

(etc) 


3-14  E12.8  Altitude 

15-26  E12.8  Axial  transition  location 

Same  as  Card  No.  2 for  increasing  altitude. 


3.1.7  Table  05  - Weather  Conditions 

This  table  inputs  the  hydrometer  environmental  conditions  for  use  in 
performing  erosion  calculations.  The  input  consists  of  certain  erosion  calcu- 
lation flags  and  a cloud  altitude  profile.  This  table  is  not  required  for 
clear  air  calculations. 


Columns 

Format 

Data 

Units 

i 

1-2 

12 

Table  No.  05 

— 

:i 

} 

1-5 

15 

NCL  - number  of  cloud  entries.  If 

| 

NCL  <0,  eititudes  are  in  meters 

i 

6-15 

E10.2 

Maximum  altitude  of  cloud 

ft  or  m 

i 

1 

1 

16-25 

E10.2 

Minimum  alaitude  of  cloud 

ft  or  m 

4 

* 

4 

t 

38-39 

12 

NOSLO  - particle  shock  layer  slow- 

— 

i 

-j 

down  flag 

0 - no  slowdown 

i 

1 - particles  are  slowed  down  as 

1 

they  impinge  on  the  shock  wave 

2 - particle  mechanical  breakup 

i 

3 

model  employed 

4 

40-41 

12 

NOHEAL  - crater  rouqhness  hsalinq 

— 

5 

j 

flag 

i 

j 

0 - no  healing  of  craters 

1 - crater  healing  by  ablation  is 

modeled 

3 

1-10 

F10.5 

Altitude  (independent  variable) . The 

ft  or  m 

3 

units  arfe  dictated  by  the  meter  flag. 

; 

!— 62) ) 

(2-62) 

(2-58) 

< 

3-16 


Table  05  (concluded) 


Card 


No. 

Columns 

Format 

Data 

Units 

11-20 

F10.5 

Particle  mass  concentration 

gm/m* 

21-30 

F10.5 

Particle  diameter 

mxl0“* 

(micron) 

31-40 

F10.  5 

Particle  specific  gravity. 

If  entered 

— 

as  zero  it  is  set  to  1.0. 

This 

table  must 

contain  at 

least  two  but  no  more  them  20  entries  and  the  alti- 

tudes 

must  be  entered  in  an 

, increasing  order. 

3.1.8 

Table  06 

- Material 

Properties 

Table  06 

is  used  to 

input  material  surface  roughness 

and  thermal 

proper- 

ties.  The  material  index  number  (MAT)  assigned  to  a given  material  need  (in 
general)  only  be  consistent  with  the  material  indices  used  in  the  nosetip  con- 
figuration input  (Table  03) . If,  however,  hydrometer  erosion  effects  are  to 
be  included  the  following  assignments  must  be  followed  due  to  built-in  values 
for  certain  of  the  erosion  correlations. 

• Carbon  phenolic  MAT  = 2 

• Tungsten  MAT  = 3 

Two  types  of  surface  roughness  are  input: 

• Laminar  or  intrinsic  roughness  (k^)  which  is  used  for  rough  wall 
transition  criteria  and  in  calculating  roughness  augmentation  to 
laminar  heating. 

• Scallop  or  turbulent  roughness  (kfc)  which  is  used  in  calculating  the 
roughness  augmentation  to  turbulent  heating. 

Turbulent  surface  roughness  may  be  input  either  as  a constant  or  calculated 
according  to  kfc  =*  K^pe~*’77  where  K^and  a maximum  roughness  height  are  input. 

Three  roughness  heating  augmentation  options  are  allowed  for;  they  are 

m No  roughness  heating  augmentation. 

• Laminar  and  turbulent  heating  augmentation  according  to  the  models 
described  in  Section  2.2.4  - but  no  hydrometer  stirring  effects. 

• Same  as  above  - but  including  hydrometer  stirring  effects. 
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Table  06  (continued) 

The  material  thermal  properties  required  include  certain  constants  (i.e., 
density,  heat  of  formation,  etc.)  plus  tabular  values  of  quantities  which  are  a 
function  of  temperature  (i.e.,  specific  heat,  conductivity,  and  emissivity) . 
Notice  that  for  the  steady  state  conduction  option  the  specific  heat  and  con- 
ductivity are  not  required  and  may  be  entered  as  dummy  values  (e.g.,  1.0). 

Card 


No. 

Columns 

Format 

Data 

Units 

1 

1-2 

12 

Table  No.  06 

— 

2 

1-2 

12 

MAT  - material  index.  If  erosion 
effects  are  included  use  the  follow- 
ing assignments. 

1 - graphite 

2 - carbon  phenolic 

3 - tungsten 

3-10 

E8 

RHO  - material  density 

lbm/ft * 

11-20 

E10.5 

TFO  - datum  temperature  for  heat  of 
"formation.  For  JANNAF  data, 

TFO  * 536“R 

°R 

21-30 

E10.5 

HFO  - heat  of  formation 

Btu/lbm 

31-40 

E10.5 

TBRPL  - laminar  blowing  rate  reduction 
parameter.  X in  Equation  (2-60) . 

— 

41-50 

E10.5 

TBRPT  - turbulent  blowing  rate  reduc- 
tion  parameter.  X in  Equation  (2-60) . 

3 

1-2 

12 

NERODE  - erosion  law  number 

— 

*].  generalized  input 

(read  the  next  card) 

*=2  carbon  phenolic  erosion 
model 

*3  tungsten  erosion  model 

Next  card  input  the  constants  of  erosion  law  (Equation  (2-62)) 
(for  NERODE  = 1 only) 


1-10 

F10.4 

A1 

— 

11-20 

F10.4 

b 

Constants  in  Equation  (2-62) 

— 

21-30 

F10. 4 

c 

“ 

31-40 

F10.4 

d 

— 

3-18 
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Table  06  (concluded) 

Columns  Format 


Card 

No. 


1-2 


3-14 


12 


E12.8 


1 


Data 


JROUGH  - roughness  heating  augmenta- 
tion  flag. 

0 - no  augmentation 

1 - roughness  augmentation,  but  no 

stirring  augmentation 

2 - roughness  and  stirring  augmenta- 

tion 

RUFL  - intrinsic  roughness  height,  k^ 


Units 


xn. 


This  table  must  contain  at  least  two  but  no  more  than  JO  entries. 
Currently  Tables  07  and  08  are  reserved  for  future  use. 


15-26 

E12.8 

RUFMAX  - turbulent  roughness  height 
flag 

in-  ) 

* 

V 

>0  - constant  turbulent  roughness 
equal  to  RUFMAX 

n 

-j 

| 

5 

• 

<0  - calculate  turbulent  roughness 
according  to 

A 

\ Klpe 

I 

fi 

t 

\ 

with  ktmax  * ABS  (RUFMAX) 

27-38 

E12.8 

K1  - turbulent  roughness  height  at 
Pe  * 1 psia.  Read  only  if  RUFMAX  < 0. 

in.  | 

t 

1-2 

12 

NC  - flag,  nominally  zero,  +1  marks 
terminal  card  of  last  material  property 
table,  -1  marks  terminal  card  of  other 
intermediate  material  property  tables. 

J 

ii 

3-10 

F10.5 

Temperature  (independent  variable) 

°R  1 

11-20 

F10.5 

Specific  Heat 

Btu/lb-°R  1 

s 

21-30 

F10.5 

Thermal  conductivity 

Btu/ft-sec-°R  1 

-*! 

0 

1 

H 

m 

F10.5 

Emissivity 

A 

3.1.9  Table  09  - Surface  Thermochemistry 

Table  09  consists  of  the  parameters  necessary  to  utilize  the  surface 
energy  balance  formulation  described  xn  Section  2.4.  The  following  paragraphs 
describe  the  input  for  a given  material. 
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Table  09  (continued) 

A single  lead  card  specifying  the  material  index  is  read  first,  followed 
by  a card  containing  control  integers,  followed  by  a set  of  tabular  thermo- 
chemistry data  cards. 

3.1. 9.1  Table  No.  and  Constants 


Col umns 


Format 


Units 


F10. 5 


Table  No.  09 

MAT  - material  index  no. 

NBPF  - flag  controlling  reading  of 
mechanical  fail  quantity  (for  use 
with  blowing  correction) 

NBPF 

0 - no  mechanical  removal 

1 - mechanical  removal  term  read 

and  used  to  reduce  blowing  rate 

CMH  - ratio  of  mass  to  heat  transfer 
coefficients  (typically  1.0). 


3. 1.9. 2 Edge  Enthalpy  Data 

Equation  (2-58)  of  Section  2.4  indicates  that  if  diffusion  coefficients 
are  not  equal  or  if  the  ratio  C^/C^  is  not  unity,  then  the  surface  energy  balance 
requires  data  about  the  edge  gases  of  the  boundary  layer.  These  data  are  pro- 
vided in  special  "edge  tables"  which  precede  each  pressure  section  of  the  sur- 
face tables  (the  various  sections  of  the  surface  tables  are  described  in  Sec- 
tion 3.1. 9. 3 below).  The  independent  variables  for  an  edge  table  are  pressure 


and 

temperature. 

Dependent  variables  are  h „ and  the  sum  T'z*  h.w. 

ew  ^ le  i 

-J 

i 

i. 

-4 

Card 

No. 

Columns 

Format 

Data 

Units 

4 

1-8 

F8.5 

Pressure 

atm 

i 

9'  16 

F8.5 

Blank 

-- 

; i 

17-24 

8X 

Blank 

— 

). 

25-33 

F9.4 

Temperature 

(*R  if  nega- 

] 

34-38 

Fi.3 

Unequal  diffusion  exponent  y 

tive  in  which 
case  enthal- 
pies below 
are  Btu/lb) 

5 
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Table  09  (continued) 
Card 


Columns 

Format 

Data 

Units 

39-47 

F9.3 

T 

Sumnation  / z?  h.w 

< re  1 

cal/gr 
(Btu/lb  if 
temperature 
is  entered 
with  minus 
sign) 

48-56 

F9.3 

Enthalpy  of  edge  gases  h_ 

Cw 

cal/gr 
(Btu/lb  if 
temperature 
is  entered 
with  minus 
sign) 

57-58 

12 

-1  (flag  signifying  that 
card  is  part  of  the  edge 
table) 

this 

gas 

59-60 

2X 

Blank 

— 

61-66 

A6 

Unused 

— 

67-78 

2A6 

Leave  Blank 

— 

5 Same  as  No.  4 for  remaining  entries  in  "edge  table"  for  this  pressure, 

(etc)  maximum  of  12  temperatures  for  each  pressure. 

Note  that  although  the  thermochemistry  codes  described  in  Section  2.4.3 
will  provide  data  decks  using  °K  and  cal/gr,  in  those  rare  cases  in  which  a 
user  wishes  to  supply  his  own  deck  and  prefers  to  work  in  ®R  and  Btu/lb,  he  may 
do  so  simply  by  introducing  a minus  sign  as  a flag  in  front  of  the  temperature 
entries . 


The  table  length  is  limited  to  5 pressure  sets  (it  may  have  only  1 pres- 
sure set)  with  not  more  than  12  nor  less  than  3 temperature  entries  in  each  set. 
The  series  of  temperature  values  may  be  different  for  the  edge  table  at  each 
pressure  set.  The  table  is  organized  as  a series  of  sections,  eacs  represent- 
ing one  pressure  and  each  preceding  the  corresponding  pressure  group  of  the  sur- 
face thermochemistry  deck  as  described  below.  The  temperature  entries  within 
ea  r ection  must  be  ordered,  either  ascending  or  descending.  Similarly,  the 
pressures  must  be  ordered  either  ascending  or  descending.  Decks  generated  by 
the  thermochemistry  programs  will  have  been  automatically  ordered  properly. 

3. 1.9. 3 Surface  Thermochemistry  Tables 

3.1. 9. 3.1  Description  of  Surface  Thermochemical  Tables 

This  table  comprises  a series  of  sections.  Each  section  represents  one 
pressure  and  one  transfer  coefficient  value.  More  than  one  transfer  coefficient 
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Table  09  (continued) 

may  be  necessary  if  the  effects  of  kinetics  on  the  surface  response  are  con- 
sidered. Nondimens ional  ablation  rate,  B^c,  forms  the  third  independent  vari- 
able within  a given  section.  The  table  has  three  dependent  variables:  SziWlliW» 
hw,  and  Tw. 

The  thermochemistry  programs  generate  separate  groups  for  each  pressure, 
one  at  a time.  All  these  groups  together  make  up  the  surface  thermochemistry 
deck.  Within  each  pressure  group  the  transfer  coefficient  values  will  be 
ordered.  Within  each  transfer  coefficient  section,  ablation  rate  entries  nec'i 
not  be  ordered  in  any  particular  way  on  the  ablation  rates;  any  necessary  order- 
ing is  made  automatically  by  the  code  as  it  reads  in  the  data. 

Users  providing  their  own  thermochemistry  decks  must  ensure  that  the 
transfer  coefficients  are  ordered,  but  the  ordering  may  be  either  ascending  or 
descending  in  each  case.  The  surface  thermochemistry  cards  are  identified  by 
a unity  flag  in  column  58,  as  described  in  the  form?4-  specification  below. 

The  number  of  pressure  groups  may  not  exceed  5 (and  may  be  only  1) ; the 
number  of  transfer  coefficient  values  in  each  pressure  group  may  not  exceed  5 
but  may  be  only  1.  If  no  kinetics  effects  are  to  be  considered  a transfer 
coefficient  of  zero  is  acceptable.  The  sequence  of  transfer  coefficient  values 
need  not  be  the  same  in  the  different  pressure  sections.  Within  each  transfer 
coefficient  section  the  number  of  ablation  rate  entries  may  not  exceed  25  e~d 
may  not  be  less  thzui  2.  The  series  of  ablation  values,  B^c  may  be  unique  or 
each  section. 

The  °R-Btu/lb  option  described  for  the  edge  tables  in  Section  3. 9. 3- 2 
may  be  used  for  these  tables  also. 


3.1.9. 3.2  Card  Formats 


Card 


i No. 

Columns 

Format 

Data 

Uni  ts 

f 

\ n 

/ 

1-8 

• F8.5 

Pressure 

atm 

t 

9-16 

F8.5 

Transfer  coefficient* 

lb/ft2sec 

: 

17-24 

F8.5 

Nondimens ional  ablation  rate 
*/PeueCM  = B’ 

— 

I 

[ 

f 

f 

f 

25-33 

F9.4 

Surface  temperature 

°K 

(*R  if  nega- 
tive in  which 
case  enthal- 
pies below 
are  Btu/lb) 

Not  provided  by  most  Aerotherm  thermochemistry  codes. 
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Table  09  (concluded) 
Card 


NO. 

Columns 

Format 

Data 

Units 

34-38 

F5.3 

Unequal  diffusion  exponent  y 

— 

39-47 

F5  3 

T 

Summation  £ziwhiW 

cal/gr 
(Btu/lb  if 
temperature 
is  entered 
with  minus 
sign) 

48-56 

F9.3 

Enthalpy  of  wall  gases  hw 

cal/gr 
(Btu/lb  if 
temperature 
is  entered 
with  minus 
sign) 

57-58 

12 

Flag  indicating  surface  thermo- 
chemistry table  entry 

— 

59-60 

2X 

Blank 

— 

61-66 

A6 

Chemical  symbol  of  surface  species. 
(ACE  and  GASKET  programs  print  such 
symbols  arranged  alphabetically  and 
truncated  from  right  end  if  necessary) 

67-78 

E12.3 

Nondimens ional  mechanical  fail 
rate  - AVp.u^  = B*ail 

— 

n+1  Same  as  Card  No.  n for  remaining  entries  in  this  section;  maximum  of 
25  entries  in  each  section. 


3. 1.9. 4 Assembled  Thermocheraical  Deck 

Figure  3-1  shows  a picture  of  an  assembled  thermochemical  data  deck 
for  several  pressures.  The  deck  corresponds  to  repeating  the  input  described 
in  Sections  3. 1.9. 2 and  3.1.9.3  for  each  pressure. 

The  surface  equilibrium  data  deck  must  be  terminated  by  two  blank  cards. 
Output  decks  of  the  thermochemistry  programs  do  not  have  such  cards,  and  the 
user  must  supply  it. 


3.1.10  End  of  Input 


The  end  of  the  input  data  deck  is  signaled  by  a single  card  with  a -1 
punch  in  columns  1 and  2. 
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Figure  3-1..  Sketch  of  surface  thermochemistry  table  make-up  for  a given 
material  including  leading  constant  cards. 


3.2 


PROGRAM  OUTPUT 


The  program  output  consists  of  output  of  the  input  data  for  check  and 
verification  purposes,  and  output  of  actual  calculations.  The  output  «f  the 
input  is  covered  in  Section  3.2.1  and  the  calculation  output  is  described  in 
Section  3.2.2. 

3.2.1  Output  of  Input 

Program  output  begins  with  an  output  of  the  restart  information.-  and 
follows  with  tne  contents  of  input  Tables  01  through  09.  Most  of  this  output 
is  fully  labeled  and  is  printed  exactly  as  input  by  the  user.  For  those  few 
output  items  not  fully  labeled , Appendix  A provides  a description. 

The  sensible  enthalpy  term  output  from  the  material  properties  informa- 
tion (Table  06)  is  defined  as 


T 

hc  » HFO  + J Cp  dT  (3-2) 

TFO 


where 

sensible  enthalpy 
heat  of  formation  (input  by  user) 
specific  neat 
temperature 

The  integration  is  performed  numerically  by  summing  over  the  table  entries. 

The  surface  thermochemistry  table  is  output  reordered  with  increasing 
ablation  rates  in  each  section.  Fcr  ec.  :h  entry  in  the  thermochemistry  tables 
the  program  computes  and  out. puts  the  quantity 


hc  = 
HFO  = 


P 

T 


TCHEM  = -h  + — 
Cw  CH 


*iw 


(3-3) 


Notice  that  this  term  combines  all  of  the  tabular  input  surface  thermochemis  ry 
terms  and  when  combined  with  the  general  surface  energy  balance  equation  (Equa- 
tion (2-59))  gives 


p U C„(H  + TCHEM)  - q . 
e e H r ^cond 


w^rad 


w w 


(3-4) 
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The  purpose  for  the  creation  of  the  TCHEM  term  is  to  reduce  computer  storage 
requirements.  For  the  simpler  case  of  equal  diffusion  and  heat  and  mass  trans- 
fer coefficients  (i.e.,  no  edge  gas  tables  and  CM  * C^)  the  TCHEM  term  is  re- 
defined to  be 

T / T T \ 

tche;;  * -hj*  + B;c  (hcw  - hww)  (3-5) 

and  Equation  (3-4)  is  still  applicable. 

3.2.2  Calculation  Output 

Two  types  of  calculation  resuits  are  output,  they  are  the  environment 
calculation  and  the  surface  energy  balance/in-depth  calculations.  The  environ- 
ment output  is  iiscussed  in  Section  3. 2. 2.1  and  the  surface  energy  balance  and 
in-depth  output  is  covered  in  Section  3. 2. 2. 2. 

3. 2. 2.1  Environment  Output 

The  environment  output  consists  of  six  tables  which  may  be  called 
according  to  the  procedure  given  in  Section  3.1.  The  following  paragraphs 
give  a brief  description  of  each  of  these  table*  and  the  sample  problem  (Sec- 
tion 4)  shows  typical  output. 

Table  1 - Summary  Information 

Table  1 contains  a summary  of  current  geometry,  trajectory  and  state 
variable  values.  This  table'  is  especially  useful  as  a quick  check  on  the  cur- 
rent nose  radius,  stagnation  point  state,  and  stagnation  point  recession. 

Also  inclu.’^d  are  the  transition  and  sonic  point  .locations,  free  stream  Mach 
number,  and  number  of  calls  to  the  environment  and  conduction  packages  (these 
are  equal  for  the  steady  state  option).  For  the  flight  option,  the  current 
value  of  altitude  is  given. 

Table  2 - Summary  Distribution 

Table  2 is  essentially  a selective  condensation  of  Tables  3,  4,  5,  and 
6.  Table  2 contains  current  body  geometry,  edge  pressure  ratio,  and  Mach  num- 
ber distributions.  Also  included  are  heat  transfer  coeff icr  •snt  distributions 
and  the  momentum  thickness  Reynolds  number  distribution.  The  LAM  flag  output 
in  Table  2 indicates  the  boundary  layer  flow  regime  and  has  the  following  three 
values : 
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• 1 - laminar  flow 

• 0 - transitional  flow 

• -1  - fully  turbulent  flow 

Table  2 finds  its  major  application  in  cases  where  it  is  desired  to  keep  the 

amount  of  printout  to  a minimum.  This  is  particularly  useful  when  a large 
number  of  output  intervals  is  expected,  and  minimum  output  is  sufficient. 

Table  3 - Entropy  Swallowing  Information 

Table  3 contains  the  current  body  geometry,  shock  shape,  and  the  entropy 
swallowing  distribution. 

Table  4 - Boundary  Conditions 

Table  4 shows  the  computed  distributions  along  the  body  of  boundary  layer 
edge  properties  and  recovery  and  wall  conditions.  Included  in  Table  4 (and  in 
Tables  5 and  6)  is  the  NTB  or  transition  flag.  This  flag  has  four  values  indi- 
cating the  flow  regime  at  a given  integration  point.  These  are: 

• -1  - laminar  flow 

• 0 - onset  of  transition 

• 1 - transitional 

• 2 - fully  turbulent  *low 

Table  5 - Heat  Transfer  and  Boundary  Layer  Quantities 

Table  5 displays  distributions  of  heat  transfer  coefficients,  heating 
parameters  and  heat  flux.  Momentum  and  displacement  thickness  and  laminar  momen- 
tum Reynolds  number  distributions  are  also  tabulated. 

Table  6 - Roughness  Heating  Quantities 

Table  6 gives  distributions  of  laminar,  turbulent  (augmented  and  smooth)  . 
and  transitional  Stanton  numbers.  Other  useful  entries  include  distributions 
of  surface  roughness  height,  both  transition  parameters  (Rek(s/6*) l//3  and 


the  turbulent  J .'ating  augmentation  parameter,  and  the  net  (laminar  and  turbu- 
lent) heating  augmentation  factor. 


3-27 


3. 2. 2. 2 Surface  Energy  Balance  and  In-Depth  Output 

Following  the  environment  output  cones  the  in-depth  and  surface  energy 
balance  output  in  that  order.  The  in-depth  output  consists  of  time  step 
information,  current  nosetip  material  configuration,  and  temperature  distribu- 
tion. The  time  step  information  gives  the  size  of  the  various  controlling  time 
steps  (see  Section  2,3.4),  the  current  conduction  iteration ,. and  problem  time. 

Since,  typically,  more  conduction  time  steps  are  required  than  environ- 
ment time  steps  the  temperature  arrays  are  output  only  for  the  last  conduc- 
tion time  step  preceding  an  environment  call.  The  time  step  information  is 
output  for  every  conduction  time  step.  For  the  steady-state  conduction  option 
no  in-depth  output  is  required  or  given. 

The  surface  energy  balance  and  new  body  point  location  output  is  the 
final  output  prior  to  the  next  environment  call.  The  output  consists  of: 

e New  body  point  locations. 

e Surface  temperature . 

e Recession  rate  (both  thermochemical  and  erosion) . 

• Nondimensional  ablation  rate  (B ' ) . 

e Crater  roughness  height. 

e Surface  heat  flux. 

e Heat  transfer  coefficient  (both  blown  and  nonblown) . 

For  the  transient  conduction  option  these  outputs  represent  the  results  of  the 
last  conduction  time  step  preceding  the  environment  call.  For  the  steady  state 
conduction  option  the  output  applies  over  the  entire  time  step  between  environ- 
ment calls. 
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Sample  Problem  No.  2 is  a steady  state  weather  flight  prediction 
of  a 9°  ATJ-S  graphite  sphere  cone  nosetip  with  a 0.65-inch  nose  radius. 

This  problem  demonstrates  the  use  of  the  flight  environment 
option.  In  addition,  the  weather  option  is  utilized.  Also,  the  sphere- 
cone  input  option  is  repeated. 
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Sample  Problem  No.  3 

Sample  Problem  No.  3 is  a steady  state  clear  air  flight  prediction 
of  a 7°  ATJ-S  graphite  sphere  cone  nosetip  with  a 0.65-inch  nose  radius. 

This  problem  repeats  the  flight  environment  option,  however,  this 
time  employing  a clear  air  condition.  Again,  the  sphere-cone  input 
option  is  used.  Also,  a short  output  option  is  demonstrated. 
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Sample  Problem  No.  4 

Sample  Problem  No.  4 is  a transient  clear  air  flight  prediction 
of  the  7°  nosetip  of  Sample  Problem  No.  3. 

This  problem  is  basically  a repeat  of  Problem  No.  3,  however, 
the  transient  in-depth  conduction  option  is  demonstrated. 
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appendix  a 


DESCRIPTION  OF  UNIABELLJD  OUTPUT  VARIABLES 


>,v  ■?’  r-" 


Name 

Description 

Units 

BPRIM 

- 

Nondimens ional  ablation  rate 

— — 

B«  - ■ 

peueCM 

CH 

- 

Blown  heat  transfer  coefficient  (i.e.  effect 
of  mass  addition  accounted  for) 

lb/ft* -sec 

CHZ 

- 

Non-blow  heat  transfer  coefficient 

lb/ft1 -sec 

CM 

- 

Mass  transfer  coefficient 

lb/ft* -sec 

DLTC 

- 

Explicit  finite-difference  stability 
limit  time  step 

sec 

DLTC  ' xi 

DLTC1 

- 

Surface  heat  flux  rise  control  time  step 

sec 

ULTC1  * “old  (^)  CTF 
\^new/ 

DLTC  2 

- 

Surface  temperature  rise  control  time  step 

sec 

DLTC2  - AtoM  §g§ 

DLTE 

- 

Time  to  next  specified  environment  call 

sec 

DLTIS 

- 

Conduction  time  step  for  first  calculation 
step 

sec 

DLTIS  = STRD  § 

DLTS 

, 

Surface  recession  control  time  step 

sec 

DLTS  = ~ — 

smax 

where  6 is  the  distance  between  the  surface  and 
adjacent  in-depth  nodes 

the 

DTH 

- 

Conduction  time  step  finally  chosen  from 
all  of  the  above  limits 

sec 

HCH 

- 

Sensible  enthalpy  of  the  solid  material 

Btu/lb 

HE 

- 

Entholpy  of  boundary  layer  edge  gases  at 
the  wall  temperature,  hgw 

Btu/lb 

HFO 

— 

Heat  of  formation 

Btu/lb 

A- 2 


r 


Name 

Description 

Units 

HZ 

- 

. _ « ^ 
Summation  2-z  xe^i 

Btu/lb 

HZW 

- 

_ * Tw 

Summation  2^z  iwhi 

Btu/lb 

RHO 

- 

Material  density 

Ib/ft* 

RIOUT 

- 

Radial  distance  from  the  body  centerline 
to  the  material  interface  boundary 

in 

ROD'.:' 

- 

Initial  surface  point  radius 

in 

SLOP 

- 

Initial  body  angle  slope  with  respect  to 
the  centerline 

deg 

STRM 

- 

Maximum  surface  temperature  rise  during 
the  previous  time  step 

•R 

TBRPL 

- 

Laminar  blowing  rate  reduction  parameter 
X sn  Equation  (2-61) 

— 

TBRPT 

- 

Turbulent  blowing  rate  reduction  parameter 
A in  Equation  (2-61) 

— 

T-DIST 

- 

Explicit  grid  temperature  array 

•R 

TT-DIST 

- 

Implicit  grid  temperature  array 

“R 

TCHEM 

— 

Ablation  parameter  defined  by  Equation 
(3-3)  or  (3-5} 

Btu/lb 

TEMP 

- 

Temperature 

•R 

TFO 

- 

Datum  temperature  for  the  heat  of  formation 

•R 

TSEN 

- 

Enthalpy  of  wall  gases,  hw 

Btu/lb 

ZOUT 

- 

Initial  surface  point  axial  coordinate 

in 

